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TECHNICAL NOTE 39^4 

AN INTEGRAL SOLUTION TO THE FLAT -PLATE LAMINAR 
BOUNDARY-LAYER FLOW EXISTING INSIDE AND AFTER EXPANSION 
WAVES AND AFTER SHOCK WAVES MOVING INTO QUIESCENT FLUID 
WITH PARTICULAR APPLICATION TO THE 
COMPLETE SHOCK-TUBE FLOW 
By Robert L. Trimpi and Nathaniel B. Cohen 


SIMMARY 


A solution to the unsteady two-dimensional laminar boundary-layer 
flow inside centered expansion waves and behind both centered expansion 
waves and shock waves is obtained by utilizing an extension of the Kdrmdn- 
Pohlhausen method. The Prandtl unsteady-boundary-layer equations are 
integrated normal to the surface bounding the flow and are transformed 
into a conical coordinate system. The resulting hyperbolic differential 
equations axe integrated in closed form for flow behind shock waves and 
by numerical methods for the flow inside or following expansion waves. 

An integral technique is applied at the discontinuities existing at the 
trailing edge of the expansion fan and at contact discontinuities (entropy 
discontinuities) so that the characteristic solution may proceed across 
these discontinuities. 

The solution to the two-dimensional unsteady laminar boundary layer 
existing at all points in an air-air shock tube is obtained by this method. 
A much shorter approximate method of solution is devised and is found to 
agree favorably with this method. This approximate method is used to 
predict the flow in hydrogen-air and helium-air shock tubes. Plots of 
wall heat-transfer rate and skin friction in air-air, helium-air, and 
hydrogen-air shock tubes are presented. 


INTRODUCTION 


Impetus to the study of time-dependent boundary layers has arisen 
because of the increased importance of the flows initiated along the 
ground and over buildings by the detonation of nuclear devices and of 
the air flow over missiles in hypersonic flight. The time-dependent 
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nature of the nuclear-shock-initiated flows is obvious; whereas hypersonic 
missile flight presents two less obvious problems, one of which is direct 
and the other, indirect. The direct problem arises because of the time- 
wise variation of the differences between conditions of the outer potential 
flow at the edge of the boundary layer and of the missile skin as the mis- 
sile encounters rapidly varying ambient conditions during its flight. To 
date, because of the relative rapidity with which the fluid boundary layer 
is able to adjust to changes, the direct problem has been treated as a 
quasi-steady one; that is, for given wall and local conditions at a time 
in a time-dependent flow, the boundary layer is equivalent to that in a 
steady flow for the same stream and wall conditions . The main apparent 
difficulty in this approach is the prediction and simulation of the cor- 
rect wall conditions since they are in turn dependent on the time history 
of the boundary layer. 

The indirect problem arises from the use of shock tubes as a means 
of experimentally simulating the very high stagnation temperatures encoun- 
tered in missile flight for purposes of obtaining data regarding heat 
transfer, skin friction, ionization, dissociation, and so forth. Correct 
interpretation of shock-tube test results requires a knowledge of the time- 
dependent flow inherent therein. The perfect-fluid flow in a conventional 
(constant-area) shock tube is well known and is depicted in figure 1. At 
time t = 0 + , the diaphragm which originally separated the high-pressure 
or driver gas (state e) from the low-pressure gas (state 00 ) is instan- 
taneously destroyed and the resultant pressure inequality is adjusted by 
the mechanism of: (l) a centered expansion wave or expansion fan (that 

is, one which originates at a single point in distance and time) progressing 
-into region e and isentropically accelerating the fluid to state £ 
following the wave; and (2) a shock wave progressing into region 00 and 
accelerating the fluid to state a behind it. (See fig. 1(c).) The 
regions £ and a, composed of the fluid originally in regions e and 00 , 
respectively, are separated theoretically by a contact surface (entropy 
discontinuity) across which pressure and velocity are constant arirl the 
temperature and density are generally different. Hypersonic simulation 
testing is usually done in region a where the stagnation temperature is 
highest. The pressure ratios across the expansion fan and shock wave are 
dependent only on the fluid in states e and °o. Detailed discussion on 
the theoretical performance of the conventional shock tube, as well as 
modifications such as increasing the area of the high-pressure chamber to 
increase the shock strength or adding a nozzle to obtain hypersonic flow, 
may be found in references 1 to 6. 

However, experiments have shown that, instead of regions cr and {; 
being regions of constant pressure and velocity as predicted by inviscid 
theory and as desired for testing, the outer inviscid flow in these regions 
is time dependent; and furthermore, it is found that the shock wave attenu- 
ates as it progresses down the tube. (See refs. 2 and 6 to 10.) These 
variations are attributable both to imperfect-gas effects and to boundary- 
layer effects. Several investigations have been undertaken to predict the 
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magnitude of the wall shear and heat transfer. In addition, attenuation 
theories based on these wall effects have been advanced in references 6, 

7, 9, and 11. 

The laminar boundary- layer problems in region a behind the shock 
have been treated in references 6, 12, 13> and 14. Reference 6 first pre- 
sented the solution to this flow by solving the boundary-layer problem, 
in a coordinate system fixed to the moving shock wave, of a semi- inf inite 
treadmill (that is, the Blasius problem with nonzero wall velocity). 
Solutions were obtained for the complete Prandtl boundary- layer equations 
on an analog computer and for the integrated momentum equation alone by 
simple computations (which, unfortunately, contained a numerical error) . 

In references 12 and 13 the equations of reference 6 were solved to a 
higher degree of accuracy on a digital computer and, in addition, velocity 
and temperature profiles were determined. The momentum equation was also 
solved in closed form by an integral method by using a fourth-order series 
for the velocity profile. Another integral solution obtained by using 
the Rayleigh velocity profile for the momentum equation and then applying 
Crocco's relation between temperature and shear was used in reference 14. 
The wall shear and heat-transfer results of these references agree very 
well. 


The situation in regard to the laminar boundary layer for the flow 
inside the expansion fan and in region £ has not been as favorable. 

Prior to the concurrent studies of the present report and the companion 
paper (ref. 15 ) > the flow in these regions was handled by rather rough 
approximations. In reference 7 > the expansion wave was assumed to be 
shrunk to zero thickness and region £ was allowed to exist from the 
leading edge of the wave to the contact discontinuity. The unsteady 
boundary layer was then assumed to be equal to that of a steady flow 
which had the same free-stream values of velocity, density, and viscosity, 
and which had traveled for the same time adjacent to a solid bounding sur- 
face. References 11 and 13 also followed the procedure of shrinking the 
fan to zero width but computed the boundary layer from the modified Blasius 
(nonzero-wall-velocity) solutions. Both solutions result in conditions 
physically unacceptable near the leading edge of the fan since they pre- 
dict infinite wall shear and heat transfer at points where the velocity 
and temperature potential approaches zero. It is also evident that serious 
errors are introduced away from the leading edge for expansion-fan pressure 
ratios not near unity because an appreciable region of varying pressure and 
velocity has been replaced by a region of constant properties. Thus a cor- 
rect solution to the boundary- layer problem outside of region cj is missing. 
The filling of this gap and the determination of the boundary layer 
throughout the shock tube by a common method was the primary purpose of 
the present report. Corner effects arising at the juncture of two walls 
or effects of opposite walls of the shock tube are neglected, and the 
boundary layer is treated as a two-dimensional unsteady flow over an 
infinite flat surface. The unqualified term f, shock tube” will be used 
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in the remainder of this report to designate a shock tube in which these 
effects are considered to be negligible. 

The theoretical investigation reported herein was conducted at the 
Langley gas dynamics laboratory. The time -dependent boundary- layer char- 
acteristic system developed in this paper for conical potential flows 
(appropriate to shock tubes and nuclear shocks) is a special case of a 
general characteristic system applicable lo the study of any time- 
dependent laminar boundary-layer flow. 


SYMBOLS 


An 


a 

B n 

b 


C n 

C w 


c 


c f 


c h 


coefficients in dimensionless velocity-profile power 
series 

speed of sound 

coefficients in dimensionless local enthalpy-profile 
power series 

ratio of thicknesses of velocity boundary layer and 
local enthalpy boundary layer in transformed 

plane, ^ 

coefficients in dimensionless total enthalpy-profile 
power series 

constant of proportionality relating absolute viscosity 

[i ^1 /*v\ ^ Tp + S 

to temperature, _ _ , / 5 — ;-§• 

ratio of thicknesses of velocity boundary layer and 
total enthalpy boundary layer in transformed 

plane, ^ 


local skin-friction coefficient, 


T w 

1 2 
2 p v u l 


local heat-transfer coefficient. 


lw 



Pw u l h l 
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C P 

c v 

_D_ 

Dt 

E 0' E lh> * * • 
E 1H’ E 2H' * • 

> 

s(l^) 

S l’ S 2h’ ' ’ * 
g 2H’ S 3H' ' ' 

H 

H* 


h 

h* 



it j 

k 

L 


coefficient of specific heat at constant pressure 

coefficient of specific heat at constant volume 

particle time derivative, + u + v 2— 

dt ox oy 

functions of | 


functions of |,t) defined in equations (E9) 


> 


functions of boundary- layer shape 


total enthalpy. 



total enthalpy difference, H - H w 


local enthalpy. 



dT 


local enthalpy difference, h - h w 



1 / 


states on opposite sides of free-stream discontinuity 

thermal conductivity 

distance normal to wall, L » 6 


l 


distance free-stream particle has traversed since 
acceleration from zero velocity by moving wave 
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M 

Wp r 


n 


0 


P 



q. 




R 

A 

R 

r 

S 

s 

T 

t 

u 


Mach number 

P C P 

Prandtl number, — — 
k 

integer denoting exponent in profile power series or 
point in solution 

arbitrary location of characteristic boundary condition 
local pressure 


heat-transfer rate normal to free-stream direction. 



dy 


local heat- transfer parameter, 



Unl 

Reynolds number, 

v w 

rate of streamwise growth of mixing region 
recovery factor 

Sutherland constant for viscosity-temperature relation 
state at shock wave 
absolute temperature 
time 

velocity component along x-axis 


u* = u - V* 
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V* 

V 



W 

X 


velocity of free-stream discontinuity in x-direction 
velocity component along y-axis 


state at wall or bounding surface 
coordinate parallel to surface bounding flow 



a 


P 


r 

7 



coordinate normal to surface bounding flow 


transformed normal coordinate for velocity profile, 

»y 

dy 


i 

A J o p e 


transformed normal coordinate for local enthalpy 

py 


profile, — 


V J 0 P 


— <iy 


local enthalpy shape parameter, B£ 

C P 

ratio of specific heats, — 

c v 


ratio of specific heats. 


value 

of 

y 

at 

edge 

of 

of 

5 H> 

6 

h' or 5 U 

value 

of 

y 

at 

edge 

of 

value 

of 

y 

at 

edge 

of 

value 

of 

y 

at 

edge 

of 


boundary layer, largest value 

total enthalpy boundary layer 
local enthalpy boundary layer 
velocity boundary layer 


slopes of positive and negative characteristics in 
|,t plane 


derivative along positive and negative character- 


istics in £,t plane, — 

w 



d-y 

d-y 

dy 


general reference state; for shock- tube case denotes 
undisturbed high-pressure region 

constant free-stream region behind expansion wave 
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A velocity shape parameter, -2A2 


P 

v 


absolute viscosity 
kinematic viscosity, 


I 

^te 


^Zs 

P 

a 


conical parameter, 

a € t 

value of £ at trailing edge of expansion wave 

value of g at entropy discontinuity 

values of g at limits of forward integration along 
characteristic lines for energy and momentum 
equations 

values of g at location of Z- shock 
density 

constant free-stream region behind shock wave 


T 


shear stress, 



A 



local skin- friction parameter, 
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X 


interval in numerical integration procedure, ln+1 " ^n 


stream function 



total enthalpy shape parameter, C 2 


GO 


transformed normal coordinate for total enthalpy 



00 


quiescent state into which shock wave advances 


1 


local free-stream state 


Subscripts not specifically designated otherwise denote evaluation 
at the appropriate state or location. 

Superscripts or subscripts + or - denote waves moving in the 
positive or negative x-direction. 

A bar under a symbol denotes nondimens ionalizat ion by appropriate 


Discussion of Validity of Prandtl Boundary-Layer Equations 

Since all the theoretical equations employed to describe the boundary- 
layer flow are based on the Prandtl boundary-layer assumptions, the valid- 
ity of these assumptions must be considered further with regard to the 
■unsteady expans ion- wave and shock-wave flows to be treated in this paper. 
The point of origin of the wave (x = t = 0) , representing, for example, 
the burst of a shock-tube diaphragm, is a singular point and must be 
excluded from the region of validity of the solution. 
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The boundary layer in the vicinity of the shock wave, since it repre- 
sents a discontinuity in the outer flow, also may be considered from a 
coordinate system fixed to the wave itself (that is, the coordinate system 
of refs. 6, 12, and 14) . The problem is then equivalent to that of the 
semi-infinite flat plate with a nonzero wall velocity in steady flow, the 
leading edge of the plate analogous to the shock in unsteady flow being 
a singularity. The Blasius-type solution for this problem is valid every- 
where except very near this singularity. Thus, in the unsteady flow prob- 
lem, the present solution, as well as those given in references 6 and 12, 
is valid everywhere except very near the shock wave. 

The ideal flow generated by centered expansion fans and shock waves 
is conical in a distance-time sense; that is, it is a function only of 
the ratio x/t. Consequently, such flows are easily handled in a |,t 


coordinate system where | = . The accompanying sketch shows the 

a e t 

significant values of <■ for the shock- tube flow of figure l(c) . 



Three regions in the expansion fan must be investigated with regard 
to the validity of the Prandtl boundary- layer assumptions. The area near 
the leading edge of the expansion fan represents a region of low velocity 
that increases -from a value of zero at the leading edge ( | = -1) • This 
situation appears analogous to that in the region near a stagnation point 
in steady flow; the Prandtl boundary- layer equations are thus assumed, as 
has been found in steady flow, to give results that are valid near the 
leading edge of the expansion fan. For an expansion fan having a finite 
ratio of leading-edge pressure to trailing-edge pressure, across the trailing 
edge of the expansion fan there- is a discontinuity in the derivatives of 
the theoretical inviscid flow. At this point the theoretical inviscid-flow 
derivatives change discontinuously from a finite value ahead of the trailing 
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edge to a zero value behind the trailing edge although the theoretical 
inviscid flow itself is continuous. An analysis of the one-dimensional 
time -dependent equations of motion at such a discontinuity shows that the 
viscosity and thermal-conductivity terms are small compared with the 
inertia and energy terms; hence, laminar diffusion is negligible. Con- 
sequently, it will be assumed that, although the Prandtl boundary- layer 
equations are invalid across such a discontinuity where the second deriva- 
tives are nonexistent or infinite, the equations are valid on each side of 
the discontinuity. A matching procedure based on the conservation of momen- 
tum and energy will then be used to bridge the trailing-edge discontinuity. 


Reference l6 has shown in regard to slip flow that the discriminating 
parameter, the square of the local Mach number divided by the local Reynolds 
number, is required to be of the order of unity or lower for the theory of 
continuous viscous flow to be valid. In the region of the fan where | 

2 

approaches , pressure, density, and temperature approach values of 

7-1 

zero ( see appendix A) , slip-flow conditions are encountered, and the 
boundary- layer assumptions are violated. The approximate values (based 
on a perfect gas, y = 1.4 and a e = 1,117 ft/sec) of this parameter for 

P+p 

various ratios of trailing-edge pressure to leading-edge pressure — ^ 

Pe 

for values of p e of 1 atmosphere and 70 atmospheres are given in the 


following table. Also shown are the values of | and of the pressure 

ratio across the shock in an air-air shock tube corresponding to the 

, -p Pte 
value of . 




Pte 

p e 

{% 

for values of of - 

1 

Shock-pressure 

ratio 

1 atmosphere 

70 atmospheres 

6.25 x 10“ 4 

8.4 x 10“ 4 

1.2 X 10“ 6 

2-9 

20 

7-5 x 10“ 6 

1.5 x 10“ 5 

2.1 X 10“ 5 

3-9 

30 

« 0 

3.4 x 10 _1 

4.9 X 10“ 5 

4.7 

4o 

~ 0 

3-4 x 10 5 

4.9 X 10 5 

5- 

44 


Mi^Vf > 

An inspection of the table shows that, based on the criteria 

U 1 


0 ( 1 ) 


for slip flow, the slip-flow regime occupies only a limited region for 
extremely strong expansion waves. 
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Outline of Derivation Procedures 

Since the development of the characteristic differential equations 
governing the unsteady-flow boundary layer is rather lengthy, a brief 
outline of the procedure followed will be presented before proceeding 
with the detailed development. The principal steps are as follows: 

(1) The Prandtl boundary- layer equations with velocity and both 
local and stagnation enthalpy as dependent variables are integrated in 
the x,y,t system from y = 0 to y = <» to obtain integral partial 
differential equations. These partial differential equations are then 
transformed to a conical coordinate system |,y,t. 

( 2 ) The wall enthalpy is assumed to be constant. 

(3) The velocity and enthalpy profiles are expressed as a power 
series in transformed normal ordinates a, p, and o>. Appropriate boundary 
conditions at the wall and free stream are applied to these series . The 

parameters and so forth are evaluated in terms of A, T, Q, b, 

c, and arbitrary constants. 

(4) Primary dependent variables Z, W, or Q are introduced and 
the derivatives of A, T, D, b, and c are evaluated. Substitution 
of these quantities into the partial differential equations results in 
pairs of characteristic equations for Z and W or Z and Q. 


Derivation of Reduced Partial Differential Equations 

The Prandtl boundary- layer equations for the two-dimensional time- 
dependent motion of a fluid over a flat plate are integrated over y in 
the x,y,t system (coordinate system fixed to the surface of the plate). 
The resulting partial differential equations obtained are then transformed 
to a conical coordinate system |,y,t. In this transformation it is 
assumed that all the outer inviscid-flow functions are functions only of 

the conical parameter | = -2—. This assumption is valid for an analysis 

a £ t 

of flows initiated by centered expansion waves or shock waves moving into 
quiescent fluid. 

The qualification was also introduced in the determination of the 
final form of the equations that the nondimensional velocity and enthalpy 
profiles could be expressed as a Pohlhausen power series in a suitable 
variable which would be chosen to account for the density variation through 
the boundary layer. 
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The following equations are the resultant reduced partial differ- 
ential equations in the conical coordinate system. (See appendixes B 
and C for details.) Because of their extreme length, they axe written 
in abbreviated form here and may be found in their entirety in appendix C 
(eqs . (Cll), (Cl4) , (C15), and (Cl6)). 


(1 . i. + 2e 2 l— — + 

V ui 9/ \ui d| 


+ 6* _t_ S9 2 _ g v w^ p w 


, u > 

2 /d — ' 


u- 


9 Ul 3t 


11 9 t 


u l\Pe / \ da / A 
' ' v 'w 


( 1 ) 


V U X 0 / Ui <U 


+ J^ + 

hj_* St 


3h v 
h x * S£ 


. + J_^ + . 


V* Jt_ S0 2 _ 2 V w/Pw\ 1_ 

T y -1 3t -l\ p ey N Pr,w\ dp 


h* 

h i* I 0 t _ jt 0 _ # 

, w V -l h l* 


( 2 ) 


_ J _§*\ +20 2 [ 2 _^ 1 + . 


36 


Ui 0 


v^i 


+ 6^; t_ S§£ = 2 

9 u-j_ St Uqlpg 



) £t 

\da , 

w 


(3) 




t 3 h w 


I 3h w 

+ — + . 


S| \ u-L >|r 


u ] _ d| %* St 


%* S£ 


□* t S\|/ 
\|r U-^ St 


2 2 —I— 2 

-1 \ p e 


1 ( 

'd - — 
H l* 

N Pr,w 

V dn> , 


St 


(4) 


w 


Equations (l) and (2) are the momentum and energy equations for u^_, hq, 
9, and 0 as dependent variables and equations (3) and (4) are for Uy, 

H^, 9, and \|r. The terms and ^ are functions only of the parame- 

ter A: — and - are functions of A, T, and b; and — and S are 
0 0 * * 
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functions of A, a, and c. It should be noted that, although A - -2A 2 
(similar to Pohlhausen), T and a are defined as equivalent to B 2 

and C 2 for the enthalpy power series. 


Assumption of Constant Wall Enthalpy 

The first step toward the simplification of this system of -equations 
is concerned with the partial derivatives of h w or H w with respect 

to | and t. If the solution is restricted to gas flows of short dura- 
tion over solid boundaries having high values of conductivity and of heat 
capacity per unit volume, it may be assumed that the wall temperature 
departs only slightly from its initial value. (See ref. 13 •) If this 
slight departure is neglected and the wall temperature is assumed to. e 
constant, all derivatives of the wall enthalpy are zero and are eliminate 
from the equations. This elimination is employed to obtain a solution and, 
furthermore, for all cases except those in which an insulated wall was con- 
sidered, it will be assumed that h w = Hw = He* 

2 

Note should be taken of the fact that, as | approaches — - m 

the fan, all derivatives with respect to | of the outer flow also 
approach zero; thus, there is the possibility in this case that the^wall 

derivative of enthalpy would be significant. As | approaches y - 

however, the boundary condition of zero fluid velocity at the wall is 
also violated as the slip-flow region is encountered and, since the equa- 
tions and boundary conditions themselves are no longer valid in this 
region, the lack of consideration of the wall enthalpy derivative is no 
considered to be significant. 


Introduction of Power-Series Profiles for 


Velocity and Enthalpy 

A transformation to a nondimensional "incompressible' ordinate . normal 
to the surface is next introduced to eliminate the density terms arising 

£. §* i V* ^ so forth . • - 


in the integrals for 


A* A’ V* V* 


A power series in 


this nondimensional ordinate is to be used to define the ^elocityand 
enthalpy profiles. Since boundary conditions applicable at infinity m 
the^phy s ic al plane are to be applied at a finite list, -- - the trans- 
formed (incompressible) plane, the assumption is made that at a finite 
distance 5 in the physical plane the boundary conditions at infinity 

This distance 6 is the conventional boundary-layer thick- 
ness of either the velocity or enthalpy boundary layers. The nondimen- 
sional ordinates for velocity, local enthalpy, and total enthalpy are a, 
p, and oo which are defined by the following relations: 
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(5) 


A = 




(6) 


V 

/ 


( 7 ) 


□ = 



In equations (5) to (7) the values of 6 are not necessarily equal 
and, in fact, for the method of solution of this paper they are allowed 
to vary arbitrarily. An alternate procedure, often employed in steady 
flow, requires the various 5 values to be equal; such a procedure appears 
to be un duly restrictive and, consequently, is rejected. Thus in the pres- 
ent analysis, the ratio of the 5 values is one of the results of the 
solution and not one of the boundary conditions enforced. Kalikhman 
(ref. 17) was one of the first to adopt this procedure for 5 U and 5ft. 

The expansion series for the various nondimens ional profiles are now 
assumed to be represented adequately by either a five- or six-term series 
as shown in the following relations : 


U 1 


u 



(8) 


— = — = Bq + B 3J3 + B 2 (3 2 + B33 5 + B4(3^ + B5P 5 

hi - h w h x * 


(9) 
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H - a 


-w 


Hi - H w 


= = Cq + Ci<jo + C 2 (u 2 + + Ci^ + C^ix? 


( 10 ) 


The boundary conditions to be satisfied for the velocity profile are 
as follows: 


Five-term series: 


As = 0 


(H) 


For a = 0, 


— = 0 
U 1 


( 12 ) 


and for a = 1, 


— = 1 
U 1 


(13) 


u 

u7 


da 


= 0 


(i4) 


d 2 -^ 

u 


— = Constant (zero except for one special case) 


(15) 


daf 
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Six-term series: 
For a = 0, 


u 

^1 


0 


and. for a = 1, 


u 



1 



da, 


0 


(17) 


( 18 ) 


(19) 


d2 JL 
^1 

cLa^ 


= 0 


( 20 ) 



da 


5 


= Constant (and later determined to be zero) 


( 21 ) 


For a = 0, 


_d 

dy 



dp _ •^ )u l 

_ _p 

dx Dt 


( 22 ) 


Similar boundary conditions sire to be satisfied for the enthalpy 
profiles with appropriate energy equations in place of equations (l6) 
and ( 22) . The reason for the assumption of a constant value (which may 
be other than zero) for the highest derivative of velocity or enthalpy 
at the outer edge of the boundary layer is discussed in the considera- 
tion of the solution at the leading edge of the expansion fan. 

The system of equations (ll) to (15) or (17) to (21) may be solved 
to yield all the A values in terms of A2 and the constant of equa- 
tion (15) or (21). (See appendix D.) The term A2 may, in turn, be 
expressed by an evaluation of equation (16) as follows: 
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But 


and 


Du^ 



S / 3u\ 

dy \ dy] 




5a U 1 


5a _ 1 

dy a p e 


(23) 


(24) 


pp 

P w^ w 


1 


(25) 


(for a given value of x, t from the assumption of the Chapman- Rube sin 
viscosity- temperature relation (ref. 18) . A combination of equations ( 23 ), 
(24), and (25) results in 



P e Pi i D&i He/Af e£ 
Pw Pw -1 ^ Pw \® j v e 


(26) 


Evaluating the left-hand side of equation (26) by differentiating 
equation (8) and evaluating the right-hand side by substituting the rela- 
tion for 5^=4 (eq. (c6)) and the definition of | yields 
Dt 


Pe Pi Pe -1 ~ ^ d -l/A\ 2 6 2 
Pw Pw Pw Ui d £ W v e t 


(27) 


A similar procedure for evaluating “the energy equations (eqs. (B4) 
and (B5)) at the wall yields: 


T = B 2 = 


N Pr h l 7-1 2 

b 2 h l* 2 1 


I 


Al - 


dpi 

dT 


d.UQ_ 

(»1 - !)Pi u l 

d| 


(28) 
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„ 1 U 1 
ft = C2 = 


c 2 2H X * 


(l - N Pr )Ai 2 + 


Nprl 


dp-j 

dT 


A 


du^ 

(up - !)p-|U 

d| _ 


(29) 


The integrals required by the definitions of 0 = / -7-II - -7-]-^- <3yj 

Jo 1\ l/ p e 

5* - f A _ dy, and so forth axe evaluated in appendix D and the 

J 0 { u iyp G 

various relations for the shape parameters are expressed in the following 
form: 


- = -(A) 

A A v 


&* 5 *m 

"a = a- ( a) 


^ = |(A,r,b) 

V V 


A , —(r) 

v v 




n* n*/ \ 


(30) 


Transformation of Equations to Reduced Hyperbolic Form 

New dependent variables , Z, W , and Q, are now introduced in 
recognition of the fact that, in laminar-diffusion problems, the dif- 
fusion distance (boundary- layer thickness) is proportional to the square 
root of the product of the kinematic viscosity and the time since dif- 
fusion began. These primary dependent variables are defined as follows: 
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z(!,t) = 


v € t 


W(i,t) . £ 


v € t 


(31a) 

(31b) 


Q(8,t) = 


v £ t 


(31c) 


Substitution of equations (3l) into equations (l) to (4) results in 
the following sets of simultaneous equations: 


) + 2z < 

S| \ Ui 0 


l 5* _i_ 

2uq 9 Uq d| 


2 + 


ai - 8 /5* _ ^1* 1^ A 

Uq y0 hq b V 0^ 


1 0 


Uq dA 


,/sa\ Jm 

Vat J " HsJ 



Ai — 
1 A 


(32) 


Sw/q _ _t *L\ + 2W 

stV -I?*/ 


q v* jl_ ^1 ^1 _ 1 dp i 

2u^ + d£ + hq* d| Pqhq d| 


1 / 1 dp l 

1 ^iNv* 

1 /- 

dpq 

(ihA 

yl\pl h l d £ 

h^* d l J $ 

»l h l*\Pl 

dT 

' de ; 


V* 

1 It 


Uq SA 




. /Sa\ 

UiJ 


-s V* 

It 


ui Sr 


4 ; 

St 1 


7Sr x 


V* r 

iiz 


»1 


Sb 


t[^ 

St, 


Hr) 

\d|/ 


1 v» /sw\ _ _2_ Vy/Pw\ 2 ^ 


+ — ~T t 
Ul <? 


St/ Uq V e \p e / V 


y - 1 u,2 “6 

Npj. b hq* 


h* 



( 33 ) 
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1^*22, 
-1 ® / 

1 s* 1 d Hl 

u l 2 U. u 5* TT * 1 □* A 

(Bi - 5) 2 + H 1 T • Hl c -o 0 

o . 

2ui 0 u-l d| 

* H ^ 



2 J 


1 0 


u-l dA 


.fe\\ JZb\ 
5t ' Ha?i 


+ — — tf — j = — — f— ) Ai — 

U;l 0 \dt J Ul V e lp e J A 


(34) 


dQ A _ J_ ^ | + 2Q < 

d| \ u-, 0 


1 

rg; , l du i , 1 

dHi 

h J_ / 

( i 

2Hl 

* ^ d| h-l* 

d| 

4 

-i 2 

\ H 1 ' 9 


Un 


2 | dp l\ 5* □ 


1_ + _L_L ^ — )C 

H x * d£ * UiHi *\ -1 d| u n 2 Pi d M A * 


/ 


Ht - 


/ 


X s o! 

2 I dp i c 0 □ + _i_ i 


U 1 H 1* u 2 p i d l At Hi dA 

Hi 

2 


dOlr 
i i 


Uj_ dft 


/oQ \ „ /dfi) 

t — - | — 

\dt J \d| j 


1 i 


Hi 


Sc 



1 O* VdQ\ _ _2_ V v/Pw\ C 1 \jr_ 
Hi ^ \dt/ Hi v e lp 6 j Npj. D 


( 35 ) 
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The partial derivatives of A, r, b, fl, and c may be eliminated 
in the following manner: Equations (27), (28), and (29) may be expressed 
in the functional form: 

A(|,t) = E 0 (£) Z(|,t) (36a) 


r(i,t) 


Eih(4) A! 2 (A) + E gh (|) A 
b 2 (|,t) 


( 36 b) 


a(i,t) 


Eih(0 a 1 2 (^) + e 2h(^) a 

c 2 (l,t) 


(36c) 


where the E terms are functions of | alone, whereas the terms Ap 

and — are functions of A only. 

A 


Since the derivatives are to be expressed in terms of the derivatives 
of the primary dependent variables, Z, W, and Q, another equation 
relating these quantities is required for each pair. The relevant equa- 
tions are the identities 


k£z_ b 2 (57a) 

W f_ 

v e t W 


z 

Q 



(37b) 
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Differentiation of equation ( 36 a) immediately gives the values for 
the partial derivatives of A: 



A_ 

Si 


1_ ^0 + 

E 0 d£ Z 8| 



( 38 a) 




( 38 b) 


and 


d ^ d log e | 

gl = 1 + 2 - — =1 + 2 A 


0 dA 
A 


d log- A 


(39) 


Differentiation of equations ( 36 b), ( 38 c), (37a), and (37 b ) and 
solution of the resultant two simultaneous equations employing the expres- 
sions above for — and — give the following relations for the deriva- 
8 t 8 ^ 

tives of r and b (or ft and c) . The equations are written for T 
and b but are also applicable for ft and c by everywhere substituting 
the equivalent functional form of the H-system for that of the h-system 
(that is, ft for T , c for b, □* for V* , Eqjj for E lh , g 2 h 

for g 2h , Z for Z, Q for W, and so forth) . 




_A 

Tb 


2 g 3h 




8 l°g e ^ 


8 log e ^ 

2 - 1 


s 2h g& V" 


^ 3 log e | 

8 log e b / Si§2h 


1 - 



8 log e b J 8 log e A 


- s x _ I ) + 2 

8 log e b / 


8 log e 1 

1 

8 log e A 

m 

E o 


dE, 


■0 


(40a) 
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r i/aw\ r 


1 /a z N 


s 2 h W W s 2 h S 1 Z \ St y 


A 

rb s 


§3h 


1 - 


N $ 

a io ge § \ a iog e £ 

X _ 2 - 1 

a iog e by a iog e a 


( 40 b) 


b l/aw\ + b 


( j ( 1 + 

2g2h W V d M 2g 2h S l Z W\ Tb 


i 

1 /az\/, . 2 A _ w ~ ee v 


a iog e 


g 3 b ” 

2 a iog e r 


+ 2 


a log e f- 

a iog e a 


t s 6h 5 l0g e V , b J 

^ a logg ^ d logg 

g 2h ^ dog e F g l g 2h 

, rb 2 g 5 h a log e r d logg 

B log e $\ ± dEg 


a log e A/Eq d| 



( 4 la) 



b l/aw', + b 

2g 2h W W 


2 g 2h giZ 


a iog e | ^ a iog e ^ 


W^Ml + 22L e Be V + o &e g 

>l z \?7 \ rb 2 S5h ^ log e r ^ log e A 


(4lb) 


where 


g. 


2 h 


i _ a iog e \ + 2 a i° 6e $ 
a iog e b a iog e r 


g 2 H 


d lo S e q ^ log e □ 

= 1 — + 2 

a log e c a log e ft 


( 42 ) 
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dAi 2 

§3h = E lh — “ + E 2h 
aA 


(43) 


s 6h 


2 ^ + A ^ 

rA H dt ■ 


(44) 


Elimination of the derivatives of the secondary dependent variables pro- 
duces the final form of the differential equations. 

For local enthalpy: 


_§_ _ U-l + 2 d lQSe A ' 
“1 ft* rr 1 d lOgp A 


5* gq' 


5 *\ 


'6z\ + L + 2 d lQSe _A t_fa ) +2 




du.q 

d| 


0 Ul " 1 / h l* 1 V* A 
2 — + II - 


5* u- 


*1 


hq b V 5*/ 


d log e A / gq\6t j 


1 + 

2 


5 * 


q d loge q d logg Eq 


Z = 


gq d log e A d l0g e | 


V w/Pw\ . 8 A 

v £ (pj ^6* 


(45) 


i «4h Vaw\ «4h>\ a.. wL/az\ + + + 

S Wi' 8 i 5h 4Wt ^ \ 4S h i* dE eA 45 /"* p A d| h i dS 

_j_A «Pj. £I.i 5 loge T d lQ8e E o . _i_ a l08e T 

hq* \pi d£ d| / A V* gq a log e A d log e | g2h a log e r 


E 0 1 Sl0g *T 

/ a log e q d log e E© 

A 

i &2h a log e r 

6h \ a log e b/ gq d log e 5 

rb 2 


\ a log e b/ 


d log 

2 L*. 2 


- d log *L 
A o e V 


d log e A a log e X 


1 3 l0 6e_J 


g 2h * l0g e b 


d log e V , 1 d l0g e E 0J A g - v - n , - 

iS ^ a loge r + gq d log e I L b 2 g 5 h a log e r d log e A a log e 


j WM 0 


v«wy vv * 


!i- . 2^_i Ul 2 da 

Npr b hq* J 0 yda / 


5 6h ' 
.2 


S log e ^ d l og e | d log e | 


( 46 ) 
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For total enthalpy: 


Ui 


e I 




1 + 2 


S* d loge A J 



— + 1+2 



d . log e A 



TT 8 * W * 10 *A 

U-l - | 2 n n q \ 


U 1 


u. 


Hn - 


C □ 0 


U l" 


/ 


0 _ 

5* 


! d l°g e Y d lo Se E 0 


g]_ d log e A d log e | 


. Vw/Pw\ 0 0 
z “ vr(p7 ) 1 A 6* 


( 47 ) 






°* «2H / \ &6 /t S 2H ' &t '| g l 


,/dz\ ./dz\ 

U /*' W 


I dpi 


+ 2 i+fc^ + ^^± + -?±Jk L_^+ J_L ^ 

I 2 ^ dg /□* Ul 2 H-j* d| H 1 *(- L d| 


Hi- — 


U 1 


2Pl d|l AD* 

Hi* 

U1 2 Pi d6 



Hi - 

— 

1 d logg Ho 

A K ( 1 _ 

a loge 

l +2 dlog ^ 

g x d log e I 

[ nc 2 6 H 

d log e C j 

' d logg A 


1 d log e d loge Ep 1 d log e ^ 


^ g 6 H 


a log, 


> e D 


d log e c 


1 SlQge T 


g2H 8 lo «e c 


8 log e i 
3 iog e n 


idlogeEoA ^ l08e p d l0g e I . 5 1°^ n] 

d log e i \q c 2 g ^ H d log e ft d. loge ^ d loge M 




( 48 ) 


where 


* n V* n V* 
O loge — ° lo Se — 

64h = 1 + 2 — — — — - “ 

o log e r o log e d 


(49) 
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8 log e 


g 


1 _ 


8 log e 

8 log e 


>. , v* 

8 log e —r 

i 4 

/ 8 log e ^ 

^ a 1 - ^ 

§2h 8 iog e r 

rb 2 5h \ 8 log e b. 

2 / 

8 log e ^ 8 log e 

0 w » „ 

V + 2 


8 log e 


i 


- 2 


- 1 


8 log e A 


2g 2h 


5 log e b 


(50) 


rb‘ 


Equations (45) and (46) and equations (47) and (48) are sets of 
hyperbolic differential equations for which real characteristics exist 
The slopes of the characteristic curves are given by the following 
expressions : 



d l°g e l 


1+2 


e 

1 5* s l " 


d log P A 

/ , , O" ^ 

/ d log e — 

1 + 2 - 

\ d log^ A 



g 4h t 

0 

^1 s 2h “ s 4h^ 


(51) 


(52a) 


for the h-method; and 



s 4H t 

-1 ^ S 2H " g 4H^ 


(52b) 


for the H-method- 
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The flow field in the |,t plane is confined to the region 


t - = S 5 < | s + = 2*1 (0 < t < ») 

a e a e 

where the waves initiating the flow have, velocities u s and u s + . 

Subsequently, it will be shown that the boundary conditions are specified 
only along i = | s “ and £ = | s + ; thus, integration along the character- 
istic lines must originate at one or the other of these values of | • 
Furthermore, as the integration progresses, even though no boundary con- 
dition is specified at t = 0 or °°, time must be increasing for the 
mathematical solution to be compatible with the physical one in that an 
event preceding another in time may influence the latter but not vice- 
versa. Thus, the requirement 


5t 



S| > 0 


determines the limiting value of % at which (in the absence of dis- 
continuities) the integration from | = £ s ” or £ = | s + must stop. If 

this limiting value of | is denoted by |*, ^ (from eqs. (51) (52a) 

o| 

or (52b)) changes sign at values of |* given by 


ij 


* = 


Si 5* S 1 


+ 2 


. , 5*' 
d log — 
c A 

d log„ A 


-1 


(55) 


= ^W* = -1 ^ s 2h s 4h 
& +* = ^Q* = Si gjf S 2H g 4H 


(54a) 


(54b) 
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These values also denote the vanishing of the coefficient of the 
term — or — in the differential equations. Singular behavior of 

the solution of a differential equation at any point where the coeffi- 
cient of the leading derivative in the differential equation vanishes 
is a common mathematical phenomenon. 


METHOD OF SOLUTION FOR THE REDUCED HYPERBOLIC 
DIFFERENTIAL EQUATIONS 


Statement of Problem 


The solution of the unsteady boundary- layer equations has been 
reduced to the problem outlined below. (The h- system alone will be dis- 
cussed, since the H-system is in general similar to it. Where major dif 
ferences between the two systems arise, these will be noted.) 


If the following five equations are given: 

(1) First differential equation for Z (eq. (45)) 

(2) Second differential equation for W (eq. (46)) 
( 5 ) A= x[l,Z(i,t)] (eq. (36a)) 

( 4 ) r = r(A,b,|) (eq. (36b)) 

(5) = §r,b,A) (eq. ( 37 a)) 


w( 6 ,t) 


w 


determine the five dependent variables Z, W, A, T, and b. The 
solutions to this set of equations with appropriate boundary conditions 
for the cases of constant free-stream flow and of a free-stream flow 
varying in a manner prescribed by a centered expansion wave are treated 
in subsequent sections. 


Method for Centered Expansion Fan 

Solution of th e problem at the leading edge .- The determination of 
the critical values of A and Z at the leading edge of the expansion 
fan (t = - 1 ) may be achieved by consideration of only two (eqs. IW 
and C4‘5)) of the five equations listed previously. The elimination of 
the other equations arises from the fact that, at h -1, the flow is 
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essentially incompressible (Mq approaches 0) ; thus, the momentum equa- 
tion is divorced from the energy equation. 


Consider first the equation for A with Eq expressed in terms 
of | by using the fact that T w = T e and = M-e : 


A = 


37-1 



(55) 


If the boundary condition is stipulated that A be finite at | — -1, 
it is evident from equation (55) that the following relations must hold at 


6 = - 1 : 


Z(-l,t) = 0 


( 56 a) 



( 56 b) 


since — (l + E) = 0 and — denotes the characteristic derivative; 
dt 5| 

that is. 





(57) 


Now the momentum equation at | = -1 reduces to 



d log e ^ 
d log e A , 



+ 


i_ —L 

u-l d| 


5* / x 

1 d loge "e"/ d lo Se E 0 [ 

gq d log e A \d log e i /£ = _i 


2A 


0 _e_ 

1 A 5* 


(58) 


52 


NACA TN 3944 


since the ratio 


h-j* 

lim — — = - (7 - 1) 
g->-l ^l h l 


Equation ( 58 ) may be reduced through further limiting and combining 
processes to 



+ 2 



== + 2 



d log e — 1 

2g e e 

g 1 d log e A I 1 + | 


37 - 1 


7 - 1 

1 - - I 


= 2Al!- 
1 A 5 * 


(59) 


or 


— 1 + 2 


d lQ6e f 1/6Z 
d log e A 



6 * 

d log e — / \ 
2 _ _ 9 (6Z ' 

g ± d log e A )_ 


2At — — 
1 A 6 * 


and finally to 




_A 

6 * 


(60) 
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Equating the values of I = -1 from equations ( 56 b) and (60) 


yields 


KI) 2 = HI) 2 £ 


or 


A = - Ai -^r 
3 8 * 


( 61 ) 


The initial value of A at | = -1 for all t must satisfy equation ( 6 l) 
regardless of the values of r and b. 

If a five-term velocity profile is assumed, substitution of the 

values of An and — in terms of A into equation ( 6 l) determines 
x A 

the equality 


120 


II 


II 

17 1 /u\ 

.. 1 

no / u \ 


A + - 

12 + — 

' 90 + 4oluJ 

9 

W) 

' 7 ct=l 


' 7 o=l 


= 0 


( 62 ) 


For A to be real, the value of (— 

l ui 


0=1 


/d 2 
U1 

dsx c 


must be such that 


0=1 


17 + 

90 4o(u ly 


1,080 


12 + ( — ’ 
u n 


> 0 


(63) 
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or 


..2 


'JL 

V u l, 


- 2l.037|^-j - 14.025 ;> o 


(63a) 


The minimum absolute value for is found (by using the equal 

w 

sign in eq. ( 63 a)) to be -0.647* A minimum value is desired since the 
velocity should approach the stream velocity asymptotically at the edge 


of the boundary layer. When this value of 


is used with equa- 


a=l 


tion (62), the starting value of A is found to be 12. 3. 

If this sequence is duplicated for the six-term velocity profile, 
the corresponding equations are: 


2kO 


7 + -(- 

6 240 \uj_ 


A + - 


10 - —f— ' 

6 \ u ly 


nt 


= 0 


(64) 


-.2 


i + -v- 

6 240 \Uj_y 


1,440 


10 - 

6\Ut I 


^ 0 


( 65 ) 


or 


It! 2 HI 

— \ + 520 /iL\ > 0 


lUn 


3 Wi 


(65a) 


The minimum value is — \ = 0, and the corresponding value for A 

u i 


is 20. 


When the energy equation was considered, the following equalities 
were employed at the leading edge : 


(a) For the characteristic slopes (see eqs. (51) and (52a)): 
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( 66 ) 


(b) For the characteristic derivatives: 


5 \ d / 5t X \ 6 _ 6 t / 6 t \ 3 

8 |J + \B|/ + \5|/_ dt 



( 67 a) 


(c) For the Z-derivatives : 


6Z\ _ 6Z 

h ) ± = as 


(67b) 


since 



is zero. 


Since 


W = 


Z 




2 


(eq. (37a)) and Z = 0, it is necessary for W = 0 to keep b finite. 
The assumption of a finite b is a physically plausible restriction 
since ii implies a finite ratio of the thermal thickness to the velocity 
thickness in the physical plane. The dominant terms in the energy equa- 
tion at | = -1 are: 
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Sta /w\ . ^5hw t *z\ + a, 


S2h \ 5 V + Si z \ a Vt 


| dh-L x d log e E q 

hi* d| g! d log e | 


■n V* 

d log e — 

0 


N V* 


6 log e 

A 

a iog e r 

rb 2 


0 


a log e — \ d log e — a log e 


6 log e ^ 

i 


i - 


+ 2 


- 2 


2 5M a log e b / d log e A a log e A 


a loge “ 7 “/ 6 log e ^ d loge — 8 loge ^ 

1 j? A V A + V 

s 2h ^ l°g e b \rb 2 S ^ h ^ log e r d log e ^ ^ log e ^ 


= — — — 

V / ^ N Pr 

(68) 


Substitution of the values already known at | = -1 and employment of 
L' Hospital's rule yields 


i • / SW \ _ 2/jA V_ 

|-,-lU/ + ~ 5W ** Np r 


(69) 


But 




lim 1 = 


- 

lim b H'if v - = b 2 /if 1(± 
VA/ Z \A y / 6Z\ 


s-»-l 


(70) 


since — = 


? 5 / + 


k) at 5 


= -1, and 


and 


lim® 

S^-lW 


l-Hil 
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Furthermore, in the limit | -* -1, equation (28) reduces to 


r 


A NPr _ A 2Np r 
2 


Combination of the preceding four equations gives 


(71) 



(72) 


This equation is identical in form to equation (6l), 


Aq A = I ^ = _ 5A2 

5* 2 


since the relations between — and A are identical when r is 


V* 


5* 


replaced by A2- 


Therefore, the initial solutions for the energy equation are identical 
to those for the momentum equation when A2 is replaced by T; that is. 


For five-term series: 


A = -2A2 = 12.3 

> 

r = B 2 = -6.i 


(73) 


For six-term series: 




A = 20 

► 

r = -10 


(7*0 


From equation (71) it is evident that the initial value of b is equal 
to the square root of the Prandtl number for all cases. 
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A like treatment of the total enthalpy method would result in initial 
values of 0 , = r and b = c. This result is expected because, at | = -1, 
Mi = 0; thus. 


Solution inside the expansion wave .- It has been proven that Z, W, 
A, P, and b are constant at | = -1 for all values of t; that is. 


of invariance with t at constant £ will be employed for the flow away 
from the leading edge. Justification for this assumption is based on the 
similarity found in the complete Prandtl boundary- layer equations derived 
in appendix E. (See also ref. 15- ) In the following paragraphs, it will 
be shown that, since the outer inviscid flow is a function only of |, 
then the boundary- layer profile form factors A, T, and b, as well as 
the primary dependent variables, Z and W, are functions only of the 
conical parameter £ . 

First, consider any unsteady two-dimensional compressible flow over 
a flat surface such that the outer potential flow is conical in the sense 


boundary- layer equations are transformed from the x,y,t system to the 
|,T],t system, the following results are obtained. (This transformation is 
carried out in appendix E.) Both the boundary- layer differential equa- 
tions (eqs . (Ell) and (E12)) and the boundary conditions (eq. (E13)) 
become explicitly independent of the time t provided that the wall 
temperature is a function of | only. Consequently, the solutions to 
the differential equations must also be independent of t in the |,T),t 
system; the parameters £ and -q are similarity parameters; and the 
solutions may be classed as similar. The velocity and static enthalpy 
profiles are then functions of only £ and q: 


lim Hp* = lim hp* 


l-»-l S-*-l 



and 



A solution having these same properties 



dimensionless parameter r| is defined as 



and the Prandtl 



( 75 ) 
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= sU,1> - <? 6 > 

he % 

A combination of these two equations shows the stagnation enthalpy profile 
to be also of the same type; that is. 


h + — . 

H(|,n,t) = 2_ = H(j, n ) (77) 

Hg H £ 

Consider the velocity profile of the integral method (eq. 8) which may 
be written in the functional form: 


u(SjH) 

UqU) 


n 



n=0 


A n (A) a* 1 


( 78 ) 


The a variable is eliminated from equation (78) by expressing, a as 
a function of n and A. (The required expression is an immediate con- 
sequence of the definitions of a and q.) Equation (78) then reduces 

to 



(78a) 


Equation (36a) is of the functional form: 


A = Eq(|) 



( 79 ) 


Substitution of equation ( 79 ) into the velocity-profile equation (78a) 
yields 


u(l>q) 

ui(S) 


n 


A n (A) 


n=0 


EnU) 


A 


-,n/2 


.n 


(80) 
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According to equation (80), the velocity profiles of the integral 
method are functions of £ , q, and A, whereas equation (75) shows the 
exact profile to be a function of only | and q. Consequently, in 
order to make the integral-method differential equations compatible with 
the exact differential equations, the boundary-layer profiles applicable 
to both sets of equations must exhibit the same similarity in respect 
to £ and t|. Thus, it is necessary that 


A = A(|,q) 

Finally, since A is the velocity form factor at an arbitrary value 
of x and t (or £ and t) and is independent of y, it is then nec- 
essary that A 4 A(q). Consequently, the following conclusion is reached: 


A( | ,q ) = A(S) 


and 


Hi) 


E 0 (l) 



E 0 (|) 



(81) 


Since all the terms of equation ( 8l) except the last one are known to 
be functions of | alone, it follows that the last term must also be a 
function of £ only; that is. 


v e t 


= Z = Z(S) 





dZ 


In an analogous manner it may be shown that W=W(|), F=r(|), 

and b = b(£) for the integral-method differential equation to exhibit 


(82) 
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the same type of similarity as the complete differential equation. 


Consequently, the solution employing 




= /dr 

dt J ^dt 


I 



0 is justified for the conical expansion-fan flow region. Simi- 


larly, it may be shown that the stagnation enthalpy variables Q, ft, 
and c are functions only of the conical coordinate | . 


Whether to employ five- or six-term power-series profiles and whether 
to use static or stagnation enthalpy methods was next determined. The 
reasoning leading to the conclusion to use the six-term profiles with 
the static enthalpy as a variable was as follows: Use of a five-term 

profile necessitates the specification of / j and 

\ Ul /Or=l 



as 


negative at the edge of the boundary layer in order to satisfy the initial 
conditions. Any profiles resulting therefrom, regardless of the variation 
of the outer flow, must be limited unduly by the restriction: 


lim — < 1.0 
owl U 1 


lim < 1.0 

ML h l* 

Consequently, the five-term profiles were eliminated for use inside the 
expansion wave. Numerical integration of any differential equation is 
difficult when the integrand becomes improper. In the integral method, 
nondimens ionalization was accomplished through division by uq_, h^*, 

or Hp*. Consequently, when these quantities become zero, the integrand 
becomes improper. In the expansion fan, these singular points occur at 
the following locations : 

At | = -1 


u^_ = 0 
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At | = -1 


hi* = 0 


and. at £ = ±1 


Hi* = 0 

The singularities at | = -1 may be easily handled analytically. 
(See previous section.) However, the singularity in Hq* at £ = 1 
would occur in the middle of a numerical or graphical procedure and would 
present a possible source of difficulty. For this reason, the static 
enthalpy system was chosen for the solution in spite of the fact that the 
stagnation enthalpy method has been more generally used in steady flow. 

Although a graphical- isocline technique could have been employed, 
the following numerical method appeared to be more advantageous in regard 

dZ 6W 

to both time and accuracy. First, two additional variables — and — 

St 8 | 

were introduced and the required two auxiliary equations were obtained 
by assuming a Taylor expansion for Z or W in the neighborhood of the 
point £ n ; that is. 


Zn±l ~ ± XZn* 



± — Zn'" + — Z n 1Y ± — Zn v + 
6 24 120 


( 83 ) 


where the plus and minus signs indicate values ahead of or behind the 
nth point, X is the interval between adjacent points, and the super- 
scripts on the Z terms express total derivatives with respect to £. 


Differentiating equation ( 83 ) with respect to £ yields 


X 2 „ in . X3 „ -i v . X 4 „ v _i_ X^ 


Z&tl = V ± *Zn" + Y Zn'" ± Y V V + ^ V ± ^ V 1 + • 


(84) 


Multiplication of equation (84) (positive signs) by — and subtracting 

the result from equation ( 83 ) (positive signs) give an equation applicable 
to the starting problem (that is, when the values are known only at one 
previous point) : 
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Subtraction of equation (83) (alternating signs) from equation (83) 
(positive signs) and addition of equation (84) (positive signs) to equa- 
tion (84) (alternating signs) yield 


Zn+l - Zn-l = 2 * z n' + ^ X^Zn'" + — X^Zn + . . . (86) 

ZA+I + ZA-l - 2 z n' + X 2 V" + X 1 ‘ z n 1V + • • • (W 


The term Z^"' may be eliminated from equations (86) and (87) to yield 


ZA+l - f(Zn + l * Z n-l) - + ^-l) + + ' • • 


1 v 4- v 


( 88 ) 


•which is applicable when the values are known at two previous points. 

Equation (88) is a form of Milne’s equation (ref. 19) but has been 
derived in a different manner to show its natural compatibility with and 
extension to equation (85). Equation (85) (and its W counterpart), the 
third derivative being neglected, was employed to evaluate the first point 
away from the leading edge of the fan and thereafter equation (88), the 
fifth derivative being neglected, was used. 

The best sequence of computational steps found for obtaining the 
solution at the n + 1 point from the nth point was determined to be: 

(1) Estimate A n+ i = A n + XA n '; A n ' obtained from equation (38a) 

(2) Solve equation (27) for Zn+q 
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(3) Solve equation (85) or (88) for = ( — j 

' n+1 

(4) Solve equation (45) for b n+1 by trial and error until the 

value of ( — ) obtained equals that obtained in step 3 

\ 5 ^n+l 

(5) Compute W n+ p from equation (37a) 

(6) Solve equation corresponding to equation (85) or (88) for 
/ 5W\ 


Wn+1 ^ 


n+1 


'5W 


(7) Solve equation (46) for W^ +1 = — , 

\B!/ 


n+1 


/ SW \ ^ 

( 8) Compare values of ( — from steps 6 and J and repeat pro- 




n+1 


cedure with new value of A n +i until the two values of 
are in agreement. 

Values of the interval X were selected as follows: 


/svA 

Wi 


n+1 


X = 0.02 (-1^5= -0.90) 

x = o.o4 (-0.88 < | 5 -0.80) 

x = 0.10 (-0.70 < £ < l.oo) 

The smaller values of X were used near £ = -1 because of the 
strong curvature of the functions A, b, and r when plotted against | . 
These functions all had infinite slopes at | = -1 fit may be shown 

that — <x (l + l) - ’ 1 ’/ 2 in the region where £ approaches -1.0 j whereas Z 

<4 

and W had finite slopes with small second derivatives. 
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Methods for Regions of Constant 
Outer Inviscid Flotf 

Simplification of the reduced hyperbolic differential equations for 
special case of constant outer flow .- The general reduced hyperbolic dif- 
ferential equations of motion and energy in the conical |,t coordinate 
system are greatly simplified for regions of constant free-stream velocity, 
pressure, density, viscosity, and so forth. 

Equations ( 36 a), ( 36 b), (37 a ), (45), and (46) represent the final 
form of the most general set of equations in |,t coordinates in terms 
of local enthalpy; corresponding stagnation enthalpy equations are ( 36 a), 
( 36 c), (37b), (47), and (48) . For the special case of constant outer 
flow these equations may be written in the following form: 

For local enthalpy: 


A = 0 


( 89 a) 


r = 


N Pr h l 7 - 1 


? 2 

Mi Ai 


Z 

W 



( 89 b) 


(90) 




2 An -1 i 

6 * A 


(91) 



+ — 1 (-) 
s 2h /W/ 1 S 2h 





+ W 





Npr 


7 - 1 2 

U-1 

b “ ^ h^_* 



da 


(92) 
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For stagnation enthalpy: 


A = 0 




2H X * 




Z 

Q 



C 2 


(93a) 

(93b) 


(94) 





Z = 2 


v w /Pw ' 


A 0 0 

Ai 

A 5* 


(95) 



f4H \ /ckA g 4H 

S 2H /\^/t S 2H W/| 


s 4h\q / dz\ 
S2H/ Z W/t 



+ Q 


^ v w /Pw\^ ^1 \|f 

vT\p 7 j *5* 5 □* 


(96) 


The following additional relations stem from the requirement that 
all free-stream quantities are constant and have been incorporated in 
the preceding equations: 


\ 

Si = 1 



s 6h " g 6H “ 0 


J 


NACA TN 3944 


47 


The momentum differential equations (91) (95) i- n both systems 

are identical and are independent of the corresponding energy differential 
equations for the case of constant outer flow. In subsequent derivations 
for this case, the equations in terms of local and stagnation enthalpy 
will be carried simultaneously. 

The slopes of the characteristic curves are: 


For equations (91) and (95) : 

t 


5* 

For equation (92) : 




t 


Si 


£_ S2h 

^ g 4h 


I 


For equation ( 96 ) : 


(98) 


(99) 



t 


* S 2H 

Uq I 

□* ®4H 


( 100 ) 


If the characteristic derivative notation of equation (57) I s used, the 
differential, equation ( 91 ) or (95) ma y he rewritten to yield 



(101) 


Since A, -, and — as well as the free- stream flow quantities are 
A 5* 

all constant, this equation may be integrated in closed form along a 
line of slope jj . After the general boundary condition that Z = 


Z 


0 
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at £ = £q is applied, the result is 


(102) 


As yet, no restrictions have been made regarding the dependence of Z 
upon t. Thus, in theory, Zq may be a function of t, and the energy 

equations may be developed accordingly. However, since the constant- 
outer-flow case is only a special case of the general conical-flow case 
discussed previously, the same arguments may be applied to show that Z, 

W, Q, A, T, fi, b, and c are functions only of £. This restriction 
immediately limits Zo, Wo* and Q q to constant values. 

It should be noted that, although the time dependency has been elimi- 
nated in direct form from the equations of motion and energy, it is still 
present implicitly in Z, W, and Q. These equations still possess char- 
acteristics and integration must proceed along the characteristic lines 
given by equations (98) to (100). 

In line with these arguments, equation (102) is differentiated 
directly; this differentiation yields 


Z - Z 0 = 



A n — — 

1 A 5* 


I 


0 


! — - 6 

1 5* 


+ Z r 


u 


0 


Hi 


*0 - S 

— - 4 

6* 


0 


dZ 




- Zo 


- = Constant 


- 1 5 ? ' 5 ° 


(103) 


After equations (89b), (90), (102), and (103) are substituted into equa- 
tion (92) and the terms are transposed, the following differential equa- 
tion for local enthalpy results: 


2 h'/'hif Al ii-Zn 
(W\ _ W v €\ p£ / 1 A 6* ^ 

W + " 2 


5? - 6 ° 


4 > g 4h 

Ui | 

®2h 




6» 


z Bl v* 


- 6 0 


Z 8 

Bl & • 5 o 


Zo - 


a^ *•*« 

v e\Pe/ 1 Lu l 


VM ^ — 

^[p e ) VV*N Pr 


i + a i~s- fti 

B 1 Al 2 J 0 J 


2 


( 104 ) 
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Similar manipulations lead to the equation in terms of stagnation enthalpy 
from equation (96): 



* g 4H , 

U-i — | 

□* 62 H 


2 ^/e»f Al a-2--zo 

Q ♦ v j }£€_1 hi * 

7 “1 n* 0 


Q 

z 


*1 


J- 

6 * 


*1 


_ 0 _ 

5* 


- 


Zq- 



2 

*w/bw\ ♦ ♦ C 1 , 
v e\pe/ GQ * N Pr 

(105) 


Equations (104) and (105) represent the alternate forms of the energy 
equations and, in general, must be solved by either a graphical- isocline 
or a numerical step-by-step procedure in conjunction with equations (102) 

and either ( 89 b) and ( 90 ) or ( 93 b) and (94). 

Analytic solution for flows generated by waves of zero thickness . - 
The boundary layer behind zero- thickness waves advancing into a fluid at 
rest permits the introduction of the relations Zq = Wq = Qq = 0 at the 

wave location (|q = Is) in the P rece<iin 6 set of equations. The familiar 

shock wave of the shock- tube flow is such a wave. Another is the mathe- 
matical model of a "negative shock" which is defined as an isentropic 
expansion wave having finite pressure ratio and zero thickness, traveling 
with the speed of sound of the undisturbed fluid, and generating a flow 
behind it otherwise identical to that generated by an expansion wave of 
finite thickness and identical isentropic pressure ratio. Hereafter in 
this report, the unqualified use of the word "shock" or "shock wave" in 
reference to a potential-flow discontinuity will designate a shock wave 
satisfying the Rankine-Hugoniot equations. Solutions to the complete 
Prandtl boundary- layer equations for the shock may be found in refer- 
ences 6 and 12 and for the negative shock in reference 13. In these 
solutions the velocity and temperature profiles are similar everywhere in 
the region of flow (that is, the profiles change shape only by a stretching 
factor normal to the wall as time increases) . The integral-solution coun- 
ter-part to such a similarity throughout the field of flow is represented by 
the constancy of A, T, b, and c everywhere in the field. 

The conditions Zq = Wq = Qq = 0 at |q = | s transform equa- 
tions (102), (104) , and (105) into the following equations: 


Z = 


2 


2 

V w/M a £ _6_ ^ ~ 

MPe/ 1 A 6* 6 . 

-1 g5f ' 


(106) 




0 ®1_ 
V V* Npr 


2 


£_ _b_ 
®1 Ax 2 



(107) 


50 


NACA TN 39 Ml 


2 A/£ 

/ Sax Q V£ V 

’«\ 2 a 0 9 

£ Al A6* 

1 

a 2 

.2 

rsM Al ii- 

2^f 

tfr* * 


if g 4H e 

z 


v «l p «/ 


5 * 


■ □* S2H 


i jl 

□ □* Npr 


(108) 


If the conditions of constancy of A, r, Q, b, and c are now 
0 5* d V* i|r r# 

imposed, — , — , — , — and ^ must be constant. It must then 

A A V V □ □ 

follow from equations (90) and (9*0 that 


Z 

W 



(109a) 


Z 

Q 



Constant 


(109b) 


Since Z is a linear function of |, the terms W and Q must also be 
linear functions of |, and their derivatives with respect to | are 
constant. Equations (107) and (108) indicate, however, that the deriv- 
atives of W and Q are functions of | unless the second terms on 
the right-hand side of the equations are identically zero. Equating 
these terms to zero in order to satisfy the similarity relationship 
produces 


l 

V 


Is V* 


Si 


+ 


A 1 N Pr 


1+2 




— — ' ) (no) 

Si A / 


ijr £s y* c^ ^1 / £ 

5 " u^ "5" + A^ Npj.1 A 


ii 

Ui A 


( 111 ) 
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Equation (110) must be solved simultaneously with those for ^ 

(eqs. (D15) and (Dl6)), ^ (eqs. (D9) and (DIO)), and P (eq. (89b)) 

to determine the values of T and b that would be required to make W 
linear. (Similar relations apply for eq. -(ill).) Elimination of the 
terms in the braces reduces equations (107) and (108) to 



W 

Z 


Vw/Pw 
v e \p 


Al il 

1 A 5* 


- U1 5 ? - 


( 112 ) 



Q 


2 


, .2 

VM 0 _e_ 

v e \P e / 1 A 6* 


Z 





(113) 


The terms W and Q then may be evaluated either directly from equa- 
tions (l09a) and (109b) or by application of the linear relation between 
these parameters and i; that is, 


w(!) 



is) 


Q(5) = (t - i s ) 


Solutions to these equations may be found for any arbitrary specifi- 
cation of free-stream and wall conditions behind any free- stream flow 
discontinuity moving into a fluid at rest. The shock and negative shock 
are special cases of this general case in that the free-stream con- 
ditions are specified by the conservation laws. 
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Solutions for the Boundary-Layer Parameters Across 
Discontinuities in the Force-Stream Flow 

Assumption of finite discontinuities to replace mixing regions . - 
In the free-stream flow generated by traveling shock or expansion waves 
moving into quiescent fluid, there are two general types of discontinuities 
which will be classified as either strong or weak* The strong discon- 
tinuity arises from a jump in any of the flow quantities such as velocity, 
temperature, density, or pressure. The shock itself and a contact surface 
(entropy discontinuity) are of this type. The weak type indicates a dis- 
continuity in the first derivatives of a free-stream flow variable. The 
flows at the leading and trailing edge of any finite expansion wave exhibit 
this type of singularity in that the derivatives change discontinuously 
from a nonzero value inside the wave to a zero value outside of it. 

In the following paragraphs, equations giving approximate matching 
solutions to the integrated conservation equations for the case of the 
weak discontinuity at the trailing edge of the expansion fan and for the 
case of the strong entropy discontinuity are presented. (The cases of 
the shock or expansion-fan leading edge are trivial.) It is to be noted 
that a line ( zero- thickness) discontinuity is employed in both cases. 
Although the entropy discontinuity actually becomes a strong entropy 
gradient either because of laminar diffusion, in which case the width 
of the gradient zone grows parabolically with time (ref. 20) , or because 
of turbulent mixing, this growth is ignored for the sake of simplicity 
since any growth that is nonlinear in time would eliminate the conical 
symmetry of the problem. The behavior of the adjusting region for the 
boundary layer at the trailing edge probably also would be parabolic in 
time, but this fact has not been proven to date. The solutions obtained 
under these zero- thickness assumptions would be only approximate but 
should give fairly good indications of the correct boundary-layer behavior 
in regions not immediately adjacent to the discontinuities. 

Matching equations for trailing edge of expansion fan .- The equa- 
tions representing the conservation of mass, momentum, and energy entering 
and leaving the discontinuity over the range of y from zero to L(L >> 5) 
are derived in appendix F (eqs. (F8), (F9), and (F10) ) . The states^ i 
and j represent states on opposite sides of the discontinuity. Since 
in this case there is only a discontinuity in the derivative, u ± = Ujj 
t). _ tj.. o = o •; and h,- = h*: and the equations take the form: 

Pi - Pp Hjj i y 


y*\ h]_* 

" ^lfe” 


V . ■* - 6-j* = | 





Vj* - Sj* 


(Ilka) 
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6i + 1 


f>>- (■- 5) - 


:-;T¥v 


(ll4b) 
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(114c) 


These equations may be rearranged to give 



(H5a) 



(H5b) 
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If i is the state inside the expansion fan, all the parameters 
of state i may be assumed to be known from the solution obtained by 
integrating from | = -1 to that point < lz* being assumed). The 

velocity of the discontinuity must coincide with the trailing edge; that 
is, v* = (u - a) i - In addition, since at j all free-stream derivatives 

are zero, Aj = 0. 


Consequently, the problem becomes one of satisfying the three simul- 
taneous equations (115) of the form = f^bj,Tj) and equation (28) in 


the form fj = P^ 


^b j) . This system is equivalent to three equations of 


e i 

— and 


u i / \ 

the form -— = f( b j ] . Such a system with only the two unknowns 
e j ' ' 

bj is redundant and solutions would not in general be expected to exist. 


In an effort to satisfy the basic equations of mass, momentum, and 
energy, another unknown variable R was introduced. The assumption was 
made that, instead of the line discontinuity, the regions i and j 
would be separated by a mixing region of thickness equal to the product 
of time and 6; £ by definition is equal to the difference between the 
leading- and trailing-edge velocities of the mixing region. In this 
manner, the conical symmetry was retained and a finite mixing region was 
permitted. In addition, storage of mass, momentum, and energy as well 
as heat-transfer and wall-shear effects was permitted inside this zone. 
Evaluation of these influences was on the basis that the values in the 
zone were assumed to be equal to the average of the values at the extrem- 
ities of the zone. The three simultaneous equations in the form 

— i = f(R,b^ were then solved. However, values of R so obtained 

9j V 0; 

located the mixing zone as extending from a point inside the expansion 
fan back to the trailing edge. Such a solution is not acceptable since 
it was assumed that Aj = 0. This solution may be interpreted as an 

indication that the physical boundary layer is influenced upstream of 
the trailing edge, since the propagation velocity for the disturbances 
is higher in the boundary layer than in the free stream. (That is, 

|_a e | > | u i — a i | - ) In other words, pressure pulses would travel faster 

in the boundary layer; thus, it is logical to expect some type of dis- 
turbance inside the boundary layer ahead of the potential-flow trailing 
edge . 


Finally, in order to obtain a usable approximation, the momentum 
and energy equations (115b) and (ll5c) were solved simultaneously whereas 
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the continuity equation ( 115a) was ignored except to provide a check on 
the magnitude of the error introduced. Such a procedure has been employed 
successfully for the momentum and continuity equations since Von Kdrmdn's 
early work on transition, in which it was assumed that at the transition 
point the laminar and turbulent momentum thicknesses were equal, conti- 
nuity being disregarded entirely (ref. 21). Comparison of the values of 

— , obtained by simultaneous solution for bi and — ^ of equations ( 115b ) 

e J 

and (115c), with that obtained from the continuity equation (using the 
solution for bj) showed a maximum discrepancy under 2 percent. Thus, 
the neglect of the continuity equation does not appear to be significant 
in this case. 

Matching equations for contact surface or entropy discontinuity . - 
The case of a strong discontinuity bounded on both sides by constant 
free-stream flow regions of the same fluid will be considered. The dis- 
continuity propagates with a velocity V* = u^ = Uj and, in addition, 

p. = p., and 7\. = A . = 0. Under these conditions, the matching equa- 

tions (eqs. (F8), (F9), and (FIO)) reduce to 

6^ = 6 j* (116) 


9i - 5 i* = 0 j ' 5 j* 


(117) 


^ ^ ^ v - 1 ^ - 


hi* hi 


hi u± 2 


2 x Ui 2 


hi UjL 2 


hi* hi , x hi hi* hi 1 

J —Ji -7 e J + -4: V - -^r 7 j' 


( 118 ) 


h 


J Ui 


Ui 


hj u-2 


Solution of the continuity and momentum equations (ll6) and (117) 
immediately gives 


0 i = 0 j 


&i* = 5j* 
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Since Ai = Aj = 0 , the following relations must also hold: 


4l = 



(H9a) 


I ± = Ii(Ai) = Ij 


(119b) 


Substitution of these equalities into the energy equation and 
rearrangement results in 



This equation of functional form f(b,r) = Constant is reduced to 
the form f(b) = Constant by use of equation (28). The resultant equa- 
tion cannot be solved in closed form but may be easily solved by trial 
and error. 


TRANSFORMATION TO AND EVALUATION OF PERTINENT PARAMETERS 


The inverse transformation from the a, p, or a) normal coordinate 
system to the physical normal coordinate y may be carried out when the 
profile functions are known. This transformation is first illustrated 
for the local enthalpy case. 

The definition of p (eq. (6)) is inverted to 


IN THE PHYSICAL FLOW 


Inversion From Transformed Coordinates to the 


Corresponding Physical Flow Coordinates 



( 121 ) 
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Substitution of the definition of h* and hi* enables equation (121) 
to be written as 


y = f p 

V p, Jo 




h, h c 


h x * 


Pe 

d3 = — 
Pi 


tV hj_* rP h * 

hT P + h c Jo h n * 


dp 


0 U 1 


(122) 


The expression for | may then be evaluated for any known distribution 

of with p. Equation (122) is suitable for transformation from p 

h]_* 

to Transformation of the velocity profile is accomplished by sub- 

V 

stituting the identity ba= p. The expression | is, in turn, related 


to a nondimens ional similar normal ordinate 




as follows: 


y _ y v _ I 

n/v V 




(123) 


Thus, 


\fz 


/hv 


\/V 


b ® p l\ h 
A X 


hq* 

— P + 

h c 


hi* r P h*_ 
h e ^ 0 h l* 


dp 


(124a) 


or 


v -k 


f ba il 


—I— ba + — / — dp I 

9 p 1 \h € he J 0 


h^ 


(124b) 
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Equations (l2Aa) and ( 124b ) are the inversion equations for the 
local enthalpy transformation from the (3- and a-plane to the physical- 
flow coordinate system. 

Corresponding equations for the total enthalpy transformations are 
as follows. (The identity ca = cjd is also employed.) 



y 

\F? 


\/z Pe 
c ®Pl 


Hw 

h e 


(JD + 




y 



\JZ Pe 

c £Pi 


H- 


w 


he 


ca + 


Hy* 

he 



Transformation of Velocity and Temperature 


Profiles to the Physical Plane 


( 125 a) 


( 125 b) 


Since the profile series result in distributions of 


u 


and 



or £ in terms of <x f or cd> it is desirable to obtain profiles 

nondimens ionali zed by an invariant reference state in terms of y, i, 
and t. In the previous section the inverse transformation of a, p, 
and cd was shown to be of the form 
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a = 



p = 



( 126 ) 


O) = 



y 



The combinations of equations (126) with equations ( 8 ), (9) > (10) 

produces the following equations for the profiles in the physical plane: 



( 127 a) 


= 1 + 






( 127 b) 


— -{y > I >b) 
He 




(127c) 


jp H pg computed by the local enthalpy method^ the following 

He 

equation with values substituted from equations ( 127 a) and ( 127 b) rather 
than equation ( 127 c) would apply: 


He 


-^y,l,t) 

be 


^e 

He 


1 

H 

JL( 

2 

V 

a e 


He 


(128) 
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Evaluation of Wall-Shear and. Heat-Transfer Coefficients 


The expressions for the wall shear and. heat transfer in basic form 
and. in various coefficient forms are derived.. In this section the designa- 
tion a after an equation number denotes the local enthalpy method and 
b, the stagnation enthalpy method. 


Consider first the wall shear: 


- ■ Hli 


_u 

Un 


T w “ Pw u l \ 


fa a) 

. da ) a=0 ' w 


Pw 1 

T w = dw u l A l ~ ^ 


Pw . 0 1 

T W - Pw “ U 1 A 1 A 




(129) 


A dimensionless shear function r £ is introduced to eliminate the 
time dependency in equations (129) with the result that t € - T e(£) and 


T w Pw Pw U 1 ^1 A 



(130) 


The rate of heat transfer from the wall to the fluid may be treated 
in a like manner: 


*w 


v w 


lb T) 

i_ 

U Pw hn*Bn * ■ 

II 

> 

\j°. 

N Pr 

P s 1 1 v j; 


.tw 


(131a) 
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1 p w -a * r 

q " = ' “ pT 1 1 D 


\/vt5 


(131b) 


^ = 




p e a e 


2\/Zl 


_ & 

1 p w Pv ^1* ^ V 

(7 - DNpr p e p e h e \/w 


( 132 a) 


^ = - 


p w p w H l* C 1 □ 


(7 - DNpr p e p e h e \[$ 


(132b) 


Analogy ratios for and f € may be simply formed and are 


A 

T _ 


± h-j* B 


bi 

Np r u-]_a € Ap 


(133a) 




N 


H x * 
Pr u l a e 


(133b) 


An alternate coefficient form which, although it is more easily 
comparable with results of others, is dependent on both |, t, and a 
flow length l may be expressed in terms of a Reynolds number based on 
flow length. The following definitions are used: 

l length potential flow has traversed since acceleration 
from zero velocity by wave 


R = 




w 


C f = 


c h = 


T W 


( 134 ) 


\ Pw u l 2 

RprQ-w 

p w u l h l* 
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Equations ( 13 I+) are combined with equations ( 1 P Q) and (131) to 
determine 




PwPw 

Pe^e 





PwPw 

P € P e 


l 

Uqt 



\|r 

Cl □ H l* /PwOw Z 
h l* /p e Pe u l t 


(135) 


( 136 a) 


( 136 b) 


The resultant analogy ratios are 


fh = ^1 

c f 2Aq 


(137a) 


fh 

c f 


cC-j_ 

2Al 


Hq* 

hq* 


(137b) 


These relations (eqs. (135) bo (137)) ar e appropriate for the 
constant free-stream flow behind a zero-thickness wave, since for this 
condition l is given by 


l 


U 1 

Uq - U S 


<X - Ugt ) 


( 138 ) 


NACA TN 3944 


63 


results and discussion 


SHOCK SOLUTIONS 


Equations (106) to (113) for the local and stagnation enthalpy 
methods were solved for the flow behind .a shock wave advancing into 
undisturbed fluid (region «) . Solutions were determined for the fol- 
lowing conditions: 

(a) Five- and six- term series represented profiles of velocity and 

enthalpy 

(b) The highest order derivative was assumed to be zero at a = 1 

and 3=1 


(c) Prandtl numbers of 0-72 and 1.0 

(d) Wall temperature equal to the undisturbed fluid temperature 

and wall temperature equal to the recovery temperature. 


Figures 2 to 4 present the results in coefficient form and show com- 
parable results for the solution of reference 12. The potential flow is 
designated by the subscript a and represents the shock-tube region 
appropriate to this case . 

The coefficients are independent of the ratio of the specific heats 
7 when plotted against a function of the ratio u s /u a . The choice of 

u s rather than — as the abscissa was made in order to represent 

u s - u a u c 

the entire range of shock strengths along a finite abscissa. 


These solutions necessitated finding the values of b which satisfied 
the equations for the local enthalpy case (eqs. ( 110) and (D15) or (Dl6)) 
the values of c which satisfied equations (ill), and the equations for 
*/□ (analogous to eqs. (D15) or (M- 6 )) when the stagnation enthalpy . 
equations were used. These equations are polynomials of high order in 
and c and, as a consequence, multiple roots were possible. Only two 
roots with real positive values of b or c were found for each case. 

One family of equations originated at the value b,c - 0 for the 
limiting weak wave and the values were much less than unity throughou_ 
the entire range: the other family originated at the value b,c - \]N pr 


for the limiting weak wave and remained of order one for all waves. The 
latter family was chosen as the correct one after application of the 
argument that, on physical grounds, an infinite ratio of velocity o 
thermal-boundary- layer thickness is not acceptable. 
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Because the momentum equation is independent of the energy equation 
when the outer inviscid flow is constant, skin friction is independent 
of Npr and is the same for both local and stagnation enthalpy repre- 
sentations. The terns Cf \[r for the five- and six-term profile series 

(fig. (2)) agree well with the corresponding values of reference 12, the 
five-term solution being ^ to b percent lower and the six-term solution 
being only from 1 to 3 percent lower than the reference values. The coef- 
ficients are significantly different from the steady- flow values of 
Cf\[R and the -unsteady- flow values are about two to three times 

as large. This difference arises because shock- initiated laminar boundary- 
layer flows are more closely related to the Rayleigh problem of the infinite 
plate set in motion impulsively than to the steady flow over semi- infinite 
plates. The flow behind the shock is equivalent to the Rayleigh flow for 
the limiting weak shock with a pressure ratio of 1.0 (ref. 12). 

Heat-transfer parameters for the special case T w - Too defined by 
equations (134) are plotted in figure 3; results of the present solutions 
are compared with results of the analysis of reference 12. As in the case 
of skin friction, heat-transfer coefficients agree well throughout the 
range of shock strengths. 

The recovery factor r is defined as 


r - ^aw ~ kl _ H av ~ hi 
’ % - h-L " - h-L 


(139) 


where h aw = H aw is the enthalpy of an insulated wall (q w = 0) . The 


A aw ~ -“aw 

condition of zero heat transfer occurs when 


v^yy 


= 0 which, in turn. 


w 


requires that or C]_ = 0. These latter will be zero when 


\ 


r,n = 6 


> 

r,n = 10 


(140) 


for the five-term series and six-term series, respectively. Although in 
figure 4 the discrepancy between the recovery factors of the local enthalpy 
methods and those of reference 12 appears to be large, the recovery- 
temperature error itself is fairly small because the Mach number is 
limited to a maximum of I. 89 . The following numerical example for air 
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(as a perfect gas) with Too = 520° illustrates the fact. At a shock- 
pressure ratio of 3 . 0 , the adiabatic wall temperature from reference 12 
is 1 443° compared with the six-term local enthalpy result of 1,411 ; 
and at a shock pressure ratio of 100, the temperature from reference 12 
is 14,737° and that of the present analysis is l4,804°. 

Insulated- wall solutions by the stagnation enthalpy method required 
recovery factors of unity for both five- and six-term profiles regardless 
of Prandtl number. Failure of this method is attributed to the low num- 
ber of terms in the profile function. In the local enthalpy method for 
five- and six-term profiles, stagnation enthalpy profiles contain 9 and 
11 terms, respectively, because of their quadratic dependence upon veloc- 
ity. Thus the stagnation enthalpy profiles may have a sufficient number 
of inflection points to yield zero heat transfer at a wall temperature 
other than that equal to the stream stagnation temperature. When five- 
and six-term stagnation enthalpy profiles axe required, however, too few 
points of inflection axe available to satisfy the boundary conditions an 
to give zero heat transfer at any wall temperature other than stream stag- 

nation temperature. The result is the unique profile ^Hoi) = 1, & H = 

anti c = o°; these values appear to be unacceptable for Prandtl numbers 
other than unity. 


Figure 5 presents velocity profiles of the five- and six-term local 
enthalpy solutions for the present analysis computed at the location 




for a value of the parameter — — equal to 2.0. Local and 

a u s - u a 

stagnation enthalpy profiles for the same case axe given in figure 6 . 

Shown also in these figures are the corresponding profiles evaluated . 
from reference 12. The integral-method profiles are slightly fuller in 
the outer part of the boundary layer but are in close agreement near the 
wall, the six-term profiles giving slightly better correlation than the 
five-term profiles. 

Solutions by the local enthalpy method with six-term profiles and 
vanishing stream derivatives were obtained for the flow behind a negative 
shock propagating into region e. Skin-friction and heat-transfer results, 
c f \jk and c h ^R, are plotted in figures 7 and 8 against | te which 1S 

an index of negative shock strength. The functions are evaiuated for 7 
eaual to 1.4 and 1.667 and Prandtl numbers of 0-72 and 1.0. As for the 
case of the shock wave, only the family of solutions with b near unity 

was considered. 


Results of the analysis of reference 13 for 7 = 1.4 and Prandtl 
number of 0-72 are also shown in figures 7 and 8(a). Good agreement with 
the corresponding case of the present solution is evident. 
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EXPANSION-FAN SOLUTION 


The initial solutions to the hyperbolic differential equations 
for Z and W indicated immediately that there was no unique mathematical 
solution to the problem; instead, an infinite number of solutions were 
possible for the six- term profiles. The reason for this behavior lay in 
the fact that the starting value of A = 20 (r = -10) not only produced 

u/h* \ 

uiW 

0 0 


the fullest profile without exceeding a value of 1.0 for 
was also the maximum point on the curves of A[ 4 I and 


but 


4 


A 1 


against A. In the range where £ approaches -1, the dominant terms 
(as has been shown previously for | = -1 in the discussion of the 
starting problem) in equations ( 45 ) and ( 46 ) contain these parameters; 
thus, it is possible to find solutions on either side of A = 20. Simi- 
larly, for each A in this region, branch solutions exist for F on 
either side of -10 corresponding to a b on each side of Np^/2. 

Furthermore, it was possible to change from one family (say A < 20, 


r > -10, b > Np r l/2J to another (a < 20, T < -10, b < Np r l/2) 


in suc- 


dT 

ceeding steps, and the derivatives obtained for — 

fairing of a smooth, although inflected, curve of r 
Thus, from a single solution at the start (point for 


and — would permit 

and b against | . 
n = 1) the possible 


solutions multiplied so that at the first step away from £ = -1 (that 
is, n = 2) there were four solutions; at the second step, eight solutions; 
or at the nth step, 2 11 possible solutions. 

The aforementioned mathematical phenomenon was eliminated on the 
basis of physical logic that such oscillating solutions were not com- 
patible (at least in the region where | approaches -1) with smooth 
monotonic variation in all the potential-flow quantities. However, 
there still remained the first four branches of the solution: 

Branch I: 


Branch II: 


A < 20; b > 


0jprj 


T > -10 


A < 20; b < F < -10 
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Branch III: 


A > 20; b > 




r > -10 


Branch IV: 


A > 20; b < y%r; T < -10 


These four branches were computed and the results of their solutions 
are presented in figures 9 bo 15- Also plotted in figures 14 and 15 are 
the results of a power-series solution to the complete Prandtl equations 
(ref. 15). This solution, obtained by expanding the velocity stream func- 
tion and a corresponding temperature-distribution function in terms of 
powers of | + 1, contains terms up to (| + l) > is exact at I 1> 
and has a progressive error as & increases. Note should be taken that 
the series solution was not completed until after the integral solution 
in the expansion fan was finished; thus, a comparison with the series 
solution was not available in the early stages of computation. 

An inspection of figures 9 to 15 reveals that, in the neighborhood 


curves of the integral solution 


of | = -1, the Z, W, T e , and q e 
are nearly coincidental in spite of the variation in the profile param- 
eters A, T, anri b. The values of t £ and q € also agree closely 

with the solution of reference 15 in this region. This near-equivalence 
of wall shear and heat transfer, concurrent with different values of A, 

T and b, is typical of integral solutions which are not concerned 
directly with the various profile shapes inside the boundary layer bu. 
instead are dependent on the integral of these profiles and their deriva- 
tive at the wall. 

As | becomes larger, branches II and IV first diverge from 
branches I and III; and the solutions to the former branches become . 
unacceptable when b becomes zero at i ~ - 0.52 since this condition 
indicates an infinite ratio of — in the physical plane. Branches I 

& u t 

and III agree with the series solution within 2 percent up to | = -0-5 
for r e and within 5 percent up to I = -0-7 for q^ 

Branch I solution does not continue past | = 0.62. At this point 
t = £ z * and further integration along the characteristic originating 
from a = -1 is prohibited. Branch III solution was halted arbitrarily 
at a » 1.0, the point at which Z approaches 0. For this solution 
| z * > 1.0 so that further integration was possible but was discontinue 
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in view of difficulties (to be discussed in a future section) encountered 
in joining the solutions to a constant-pressure region. 

Velocity and enthalpy profiles of the four branch solutions and the 
power-series solution are shown in figures l6 and YJ for values of 
| = -0.9, -0.7, and -0.5 (no solution for branches II and IV at | = -0.5)* 
All solutions agree very well in regard to velocity profile (the over- 
shoot in velocity in branch III at these values of £ is almost negli- 
gible) and the agreement at <■ approaching -1 between the integral and 
the exact (that is, the power series as | approaches -1) is much better 
than that for the shock solutions. (Compare figs. 5 and l6.) Although 
the enthalpy profiles of branches I and III agree closely with each other 
and the series solution, those of branches II and IV indicate a different 
and very improbable behavior from the beginning. The profiles for 
branches II and IV show minimum values of h (maximum values of h*/hi*) 

to exist inside the boundary layer, and such a minimum does not appear to 
be physically logical since the convective derivative of free-stream tem- 
perature is always negative whereas the wall remains at a fixed temperature. 

Of the two remaining branches, I and II, arguments could be advanced 
originally for the retention of both. For the case of A < 20 (branch I) 

Dgl 

consider the argument: Since the convective pressure gradient 

D log e t 

27 

increases monotonically from a value of - a-t | = -1 to a value 

7+1 

of 0 at | = — - — (value of | for p x = o'), from steady-flow 
7 - 1 ' ' 

analogy A would be expected to decrease monotonically. For the case of 
A > 20 (branch III), the following reasoning can be applied: 

(a) If the solution were to be continued for an expansion fan of 
infinite pressure ratio ^p-^(l) = 0, | = — , consider the expres- 

sion relating A with Z: 


A = 



1 

1 + I 


37-1 


7+1 
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For A to be finite at | = 


7-1 

but also that all n derivatives 


requires not only that Z is zero 
r\ n 7 - 1 

for n < — — 1 are zero. 

n - 7 - 1 


= 0 to be true for 


d! 1 

For A to be zero would also require 

d n z 

These restrictions on — — are very rigid and a violation would not be 

d| n 


d^Z 

d| n 


n<52L^J:. 

7-1 


unexpected $ thus, A = 00 ^when | = 


is possible. 


(b) From steady-flow experience, when a strong favorable pressure 
gradient is imposed on a flow in which the wall is hotter than the free 
stream, the local free-stream velocity is often exceeded inside the 
boundary layer (that is, A > 20) . 

The subsequent completion of the series solution indicated a pref- 
erence for the case A < 20 since in the vicinity of i = -1 the 

profile — never exceeds 1. (See fig. l6.) It should be noted, how- 
U 1 

ever, that, in the region | > -0.5, the profile of reference 15 shows 
a tendency to exceed 1.0. This value of i is such that the error 
arising in the series solution by neglecting 0(1 + |)^ is of the same 
order as the velocity overshoot. 


In view of the foregoing arguments, the following policy was adopted: 
The branch of the solution for A S 20 is considered to be the correct 

and most accurate one. Results of this solution including pertinent 
derivatives may be found in table I. Results for both branches I and III 
will be presented in some cases of t £ and qg since, as mentioned 

before, the integral method itself cannot distinguish the behavior inside 
the boundary layer. 


From an inspection of figure 15 , a maximum value of is evident 

in both the integral and power-series solutions. Such a maximum was 
expected since the ratio of free-stream total enthalpy to wall enthalpy 

2 

| = -1, decreases to 


has an initial value of 1.0 at 


at 


7+1 

and increases thereafter so that it passes through 1.0 again at 
(See accompanying sketch.) 


I 

I 


0 , 

1. 
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Since the difference in these enthalpies may be considered a measure of 
the influence of the local free- stream flow on heat transfer, the heat 
transfer to the fluid would be expected to increase from an initial value 
of zero to a positive maximum, return to zero, and finally become nega- 
tive. The fact that the positive maximum at £ = -0.6 and the return 
to zero at i = 0.5 for solution I occurred at lower values of £ than 
the maximum and zero difference, respectively, in the stagnation enthal- 
pies reflects the past history of the boundary- layer flow. Energy has 
been added to and stored in the boundary layer so that at a given value 
of i the total enthalpy distribution in the boundary layer is entirely 
different from that which would have existed if the local free- stream 
condition had been constant to that point. Similar behavior occurs for 
branch III. 


A maximum also exists in the curve of the variation of t € with f- 

for the integral solutions (fig. lb), (it does not appear in the series 
solution since this effect requires consideration of terms of 0(1 + 
or higher.) A simple explanation for its occurrence in this particular 
range of | is not obvious although the existence of a maximum at some 
point in the expansion wave will be argued. 


Although the computations have not been carried out to large values 
of | (and, as discussed previously, the no-slip conditions are violated), 
an idea of the behavior of r e and under the Prandtl boundary-layer 

assumptions can be obtained from physical considerations. For an infinite 
expansion-fan strength = 0^, the following conditions apply: 

2 

*1= T^T 


% 

Hw 



0 + 


2 


7 - 1 


NACA TN 3944 


71 


Thus, although the free-stream density, pressure, and temperature 
have 1 become zero, the nondimens ional free-stream velocity u-j_ and 

H 1 

stagnation enthalpy ratios — are finite. If this situation is repre- 

H w 

sented as an infinitely thick boundary layer separating a surface from 
a stream of finite velocity and total enthalpy, the wall shear and heat 
transfer would then become zero. The following sketch indicates the 
anticipated form of the curves for r e and cjg, the solid lines denoting 

computations and the dashed lines, possible extrapolations. 



2 

For the wall shear to approach zero as | approaches ■ _ -£ requires 

a maximum on the curve of the variation of t £ with g. Similarly, the 
curve for the variation of q^ with | must have both a maximum and 
minimum value. 


Dimensionless velocity and enthalpy profiles inside the expansion 
fan for various values of | are shown in figures 18 and 19 for "the 
branch I solution. One may note in figure 18 that the value of the 
dimensionless normal ordinate at the edge of the velocity boundary layer 

y increases nearly linearly with g in the range -0.4 S g ^ 0.6. 


_Jk - 2 ( i + |) and A does not change drastically, it is to 

a e 7+1 



Since 
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be expected that the slope of the velocity profile at the wall would be 
nearly constant in this range and hence result in an approximately con- 
stant skin-friction parameter. (See fig. 14.) 


Although the value of 



is also approximately linear in 


the same range, the wall heat transfer is not constant because: 


(a) 


h}* 

he 


is not linear in £, and (b) r approaches and passes through the value 10. 
Condition (b) requires an inflected enthalpy profile and a resultant 
reversal of the sign of the heat transfer. (See curve for | = 0.6 in 

fig. 19.) 


The behavior of the solutions at the leading edge of the expansion 
fan merits consideration. Consider the following limiting equations: 

For six-term series: 


lim Z = lim a /— \ (| + 1) 

1—1 1—1 w 


lim 0.1148(1 + |) 

1—1 


(141) 


For six-term series: 


^ Up An U*i 

lim r £ = lim — — == ^ = lim 1.118 — •! — ■ ■ (l42a) 

£_>_! i— 1 \Ji + 1 £-»-! \Ji + | 

It 

For five-term series, [ — ) = -0.647: 

\ Ul /cv=l 


lim Tg = 

£— 1 


Si 



Si 

lim 1.125 — 

i^-i \Ji + i 


lim 

s — 1 


(142b) 
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Consider the same limiting equations for the case of a vanishingly 
weak positive or negative shock (that is, us -> t ae, Ui -» 0) : 

For the six-term series: 

Hm Z = 2A-L f X lim (±| s + i) = 0.20406 lim (±| s + l) 
t s -»±l A Is""* 1 ^ ±IL 

±1 |-* ±1 !-» 


For the six-term series: 


lim 

Is - * ±1 
!-» ±1 


r e = lim 


Is - * 



l-» ±1 


— = lim 


0 . 339^1 


Is - * 11 \/±| a + I 

±1 v 


( 145 ) 


( l44a) 


For the five-term series: 



lim 

! s ^±l 

l-> ±1 


0 . 547^1 


(l44b) 


The ratios of the limiting values in the expansion fan to those in 
the shocks are then 

For the six-term series: 


^exp 

lim 

|-»-l Z s hock 


0.1148 

0.204o6 


O .5626 


(145) 
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For the six- term series: 


(?.) 

\ fc /exp 

rn — 

\ e / shock 


1.118 

0.559 


2.000 


For the five-term series: 


lim 

4-»-l 


( Tg )exp 

ttt 

/shock 


1.125 

0.547 


2.06 


( l46a) 


(l46b) 


Similar ratios may be found for W and q € . The power-series 
solution, which is exact at % = -1, also yielded a ratio of r £ exactly 

twice the Rayleigh (vanishingly weak shock case) value. These results 
will be applied in a later section to obtain approximations to the inte- 
gral solution for any gas. 


DISCONTINUITY SOLUTIONS 


At Trailing Edge of Expansion Fan 


Solutions to equations (ll4b) and (ll4c) for branches I and III 
obtained at the weak discontinuity arising at the trailing edge of the 
expansion, fan are shown in figures 20 and 21. For type I solutions, the 

ratio — (the subscript i indicates conditions inside fan at trailing 

Q j 

edge; the subscript j, conditions outside fan at trailing edge) was in 

0 • • 

the range 1.05 = — — ^ 1.12 whereas -r- was in the interval 

9 j Pj 

O.93 < < 1.10. For branch III, the ranges were 0.56 = ~ = 1.12 




and 0 98 A — 0 2.2, respectively. As mentioned previously, the values 

' /j ‘ 


from the continuity equation (113) varied from those given in figure 20 
by no more than 2 percent. For branch III and for values of > 0.6, 

no solutions could be found with a physically acceptable value of b 
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(that is, 0(1 - 10)). This inability to match the values at the trailing 

edge may be construed as a further argument against the use of branch III 
solutions. 


Since 





the discontinuous change in 


Z 


and W at the trailing edge represents a Z-jump and a W-jump in the 
boundary- layer characteristic system. These jumps may be considered 
analogous to the change in values of the Prandtl-Busemann parameters 
(Ci and C2 of ref. 22) at a vortex sheet in the steady-flow character- 
istic system since there is a discontinuous variation without any inter- 
section of characteristics of the same family. Of course, these Z- or 
W- jumps exist only mathematically, since in the physical flow viscous 
effects would diffuse the jumps into steep gradients with the peaks 
rounded off, a situation similar to that in the steady flow where the 
vortex sheet would undergo viscous diffusion. 


(^e)i (3e)i 

In figure 21 are shown the ratios — — and — of the shear 

l T 0j («.j 

stress and heat transfer on opposite sides of the trailing-edge discon- 

A 

tinuity. Both solutions exhibit p* ^1.0, a relation which was expected 

T j 

for a flow emerging from a region of favorable pressure gradient into 
one of zero gradient. 


For weak expansion waves the trend of the heat transfer at the 
trailing-edge discontinuity is that which would be expected from Reynold s 
analogy in that the heat transfer decreases similarly to the skin fric- 
tion. However, for values of > -0.2, the opposite occurs and the 

heat transfer in region £ at the trailing edge exceeds that inside the 
expansion fan. No obvious argument is apparent to justify this behavior. 

Velocity and enthalpy profiles on opposite sides of the trailing- 
edge discontinuity are shown in figures 22 and 23« Corresponding pairs 
of velocity and enthalpy profiles show only a slight deviation across 
the discontinuity. 

The velocity profiles inside the fan are fuller near the wall than 
in the region £ because, as mentioned previously, the favorable pres- 

9 • 

sure gradient is absent in £. Since -i « 1.0 and the thicknesses of 

9 j 

comparable pairs of velocity profiles are nearly equal, the adjustment 
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of the boundary layer at the trailing edge of the expansion fan is 
essentially a transfer of momentum from the inner part to the outer part 
of the boundary layer. 


Equation (120) was solved for various expansion- and shock- wave 
pressure ratios for the special case arising when the same gas is used 
in both the high-pressure and low-pressure chambers of a shock tube. 

These solutions indicated no change in shear stress across the discon- 
tinuity, and only a small decrease in the magnitude of the heat-transfer 
rate although W, F, and b show extreme discontinuities which, of 
course, would be rounded off by viscosity. At first glance this heat- 
transfer result might appear to be incompatible with the physical flow 
since across the contact surface there is always a discontinuity in the 
free-stream stagnation enthalpy. For a recovery factor near 1.0, the 
heat-transfer rate might be expected to vary as the difference between 
the free-stream stagnation enthalpy and the wall enthalpy. (For example, 
when = 1.0, this difference is -O.rThg in region £ and O^hg in 

region 0 ) . However, on further inspection, it appears physically plausi- 
ble to expect the heat transfer to change only slightly although the free 
stream changes markedly. As the outer contact surface progresses at free- 
stream velocity, the colder fluid forming the outer part of the boundary 
layer just behind the contact surface overrides the hot fluid near the 
wall. Consequently, since the fluid adjacent to the wall originally came 
from region a with a high stagnation enthalpy, it still possesses enough 
energy to maintain a large rate of heat transfer simultaneously to the 
wall and to the outer part of the boundary layer. The enthalpy and veloc- 
ity profiles for the case = 1.0, |-(- e = 0, = 0.833 in "the region 

near the entropy discontinuity are shown in figures 24 to 26. These fig- 
ures show the initial hot-fluid entrainment near the wall and the subse- 
quent rapid readjustment as the energy of this fluid is given up to the 
wall and outer boundary layer until finally all identity with the previous 
flow in region 0 is lost and the boundary layer behaves as one would 
expect if region £ had been the only influencing factor. Although the 


mately equal on both sides of the discontinuity (see fig. 24), the thermal 
boundary-layer thickness immediately behind the discontinuity is many 


At the Contact Surface in a Shock Tube 


velocity boundary- layer thicknesses / at 




are approxi- 


times that ahead of the discontinuity 



= 0(10)\ . A 
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much larger thickness of the thermal layer in this region is justifiable 
on the physical basis that the temperature is forced to go first from T w 
to a value near T a > T w in the hot sublayer and then revert back to 
Tf < T w in the outer layer. (See figs. 25 and 26.) No such behavior 


is required for the velocity since U£ * u 0 . Of course, because of the 

effect of mixing and diffusion, the discontinuity of zero thickness itself 
is a physical impossibility; thus, the profiles in the immediate neighbor- 
hood of the contact surface do not physically exist, although the trends 
exhibited are probably correct. 


Since the unifying influence on all the various particular solutions 
to the localized phenomena existing in the shock-tube flow (the effect of 
only one wall being considered) is the characteristic system, brief men- 
tion of a few significant factors will be made at this point. First, all 
solutions have been obtained in conical form in that everywhere Z, W, 

A, T, and b are functions only of |. Consequently, although the 
slopes of any characteristic curve are equal to the product of time t 
and a function of |, the limiting lines dividing the integrations from 
| = -1 and those from I = l s are invariant with time. (See eqs. (53) 

and (5ka).) A boundary- layer characteristic diagram for tte = 0 is 

shown in figure 27 and may be considered typical of the diagrams for all 
! te S 0.5. Free-stream flow characteristics and a particle path are also 

shown. The regions -1 £ S £ ! z * and | w * ^ I ^ ! s emerge as two com- 
pletely independent regions since no boundary- layer characteristic in 
one region was generated in the other. Only the region 

depends on both the expansion- and shock-wave strengths because char- 
acteristics which originated at | = -1 and | = | s cross this region. 

From the characteristic theory, in the region of validity of the 
equations, there can be no discontinuity in the functions Z or W but 
there may be discontinuities in the derivatives of Z and W. At the 


DISCUSSION OF CHARACTERISTICS SYSTEM 


The Characteristic Plot for the Shock Tube 


location £ = | z *, there are discontinuities in both 



and 




whereas at the location £ = only 


is discontinuous. 
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Location of | w * by the "Anchor Point" Method 

The fortunate circumstance that Z is continuous across the entropy 
discontinuity (for the case of the same gas on both sides of the discon- 
tinuity) permits the analytical determination of both the location of 
| w * and values of W, T, and b along | w * (all of which would other- 
wise have to be obtained by numerical and graphical methods) and thus an 
"anchor point" is available toward which the graphical- isocline methods 
used for | t < | £ | d must converge at | w *. The mathematical device 

employed to obtain this result is to shrink region o to zero, to let 
region £ exist from |^. e to | s , and then to solve the resulting 

positive-shock equations. This distortion of the potential flow leaves Z 
completely unchanged from its distribution when a is nonzero, since Z is 
still linear and continuous in the interval | ^ Is • ( See accom- 

panying sketch.) 



Furthermore, the curve for the variation of W with | will remain 
unchanged in the interval | ^ | w * since the W- characteristics orig- 
inate at I = -1 and the Z-characteristics, although they originate at 
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| = | s , carry the unchanged (correct) value of Z at each £. In the 
interval | > |y*, the distribution of W is now both linear and con- 
tinuous but does not represent the true distribution) however, at | = £ w *> 

the value of W is correct and can be computed analytically from its 
linear variation from 0 at | = The terms P and b, which are 

constant for | w * ^ | S | s in the distorted flow when region a has zero 

width, are also the undisturbed values at £ = £y*. 

Physical Significance of the Critical Characteristics 

The mathematical singularities 321(1 ^w* of charac ‘ teris " tic 

integral solution have been found to be the limits of forward integration 
along a characteristic if the mathematical sequence of events is to follow 
the physical sequence. The further significance of these integral singu- 
larities in relation to singularities arising in the complete Prandtl 
boundary- layer equations and in the physical flow will be discussed in 
this section. 

The complete unsteady Prandtl boundary- layer momentum and energy 
equations given in appendix E (eqs. (Ell) and (E12) ) in the transformed 

coordinate system .£,T),t contain the term £ - as a coefficient for 

d f and g. Consequently, singular behavior might be expected 

dq d£ d£ 

, df 

in the region where £ approaches — . 


location in the boundary layer at time t of particles which were orig- 
inally at x = 0 and t = 0. Now at x = 0, t = 0, when the flow was 
originated by the instantaneous formation of the expansion wave or shock 
wave or both, a singularity was propagated noamal to the wall out to all 
values of y. This singularity coincided with the expansion wave (which 
at t = 0 has zero thickness regardless of strength) or shock wave or 
both. In the case of the shock tube when both expansion waves and. shock 
waves appear simultaneously with the diaphragm burst, the singularity^ is 
also concurrent with the position of the contact discontinuity at t = 0 

The boundary-layer assumptions admit no diffusion in the streamwise 
direction so that this singularity is conveyed along with the boundary- 
layer flow and, at any time t, will appear at a value of 


Such singular points 



are roughly related to the 
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x 



u (x,y,t) 


dt 


Consequently, in the complete boundary- layer equations, the mathematical 

singularity at £ = — may be assumed to have a correlation with the 

dT] 

physical location in the boundary layer of the generating singularity or 
contact discontinuity. In the shock-tube case the singularity reflects 
the separation within the boundary layer of the fluid accelerated by the 
shock from that set in motion by the expansion wave. The contact surface 
is analogous to the "critical characteristic" for the potential flow. 

The y- location of this singularity will vary from y = 0 at x = 0 to 
y = 5 at x = ujt. 

Since integral solutions are not affected by details of the flow 
inside the boundary layer but only take into account integrated average 
effects, the aforementioned singularity will also appear under averaged 
conditions. Instead of the singularity distribution, at various values 
of y for various values of £, which is found in the physical flow and 
in the complete Prandtl boundary- layer equations, the singularities are 
lumped together and appear for all values of y <: 6 at a particular 

£ = The fact that the averaged location of the singularities in 

velocity and enthalpy do not coincide in the integral solution (that is, 
£ Z * 4 £ W *) ma y be attributed to the fact that the velocity and enthalpy 

boundary layers are identical neither in profile shape nor in thickness. 
Thus the two critical characteristics represent different averaged posi- 
tions of the same physical singularity. 

Additional discussion of the appearance of singularities in the com- 
plete boundary- layer equations for the expansion fan may be found in ref- 
erence 15. In reference 23 the case of a semi-infinite flat plate impul- 
sively set in motion is treated, and the analogous problem of the relation 
between a similar singularity arising in the complete equations and that 
in the integral equations is examined. 


Momentum and Energy Shocks 

In any characteristics system, the merging of characteristic curves 
of the same family denotes the end of the continuous distribution of the 
dependent variables. For potential flow, the resultant discontinuous 
change is called a shock wave and this definition will also be applied 
to any discontinuous change in the boundary-layer parameters Z and W 
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resulting from the intersection of characteristics of the same family. 

Such a shock is formed by the Z- characteristics for expansion waves with 
values of He ^ 0*508 because the slope of the characteristics orig- 
inating from | = -1 is greater than the slope of those originating 
from £ = | s . In figure 28 characteristic diagrams for the Z- characteristic 

family are sketched for the cases of 


(a) S te = “°-5 « 0.508 

(b) | te = 0.4 -> 0.508 

(c) He = 0-508 

(d) O.508 < | te ^ 0.62 

He = 0.55 

(e) and (f) He = ^-00 


The relative slopes of these sketches are distorted in order to indicate 
more clearly the convergence or divergence or both convergence and diver- 
gence of the characteristics. For values of He < O.508 (fig. 28(a)), 
there is no merging of the characteristics, although as He “* O.508 
(fig. 28(b)) the heavy density of characteristic lines behind the 
expansion-fan trailing edge indicates a strong gradient. When | t e = °-508 
(fig. 28 (c)), the Z boundary- layer characteristic family in regions £ 
and a has a slope equal to that of the upstream potential-flow char- 
acteristic (expansion-fan trailing edge) and all the Z-characteristic 
lines from | = -1 are trapped at He = °-5° 8 311(1 extend to t = » 

along He = O.508. This condition is the initial formation of the 

Z-shock due to the merging of the characteristics. For values of He 
between O.508 and 0.62, the sketch (fig. 28(d)) is drawn for convenience 
with the Z-shock located concurrent with the expansion-fan trailing edge, 
but this location has not been established mathematically. The only math- 
ematical restriction on the Z-shock location under this condition is that 
the slope of the shock must be between the slopes of the bounding 
characteristics. 
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Methods for Extending Present Results 

Since \ ^ 0 in the fan, it can be proven for £ te > 0.508 that 

any Z- characteristic entering the fan must be bent from its direction in 
region £ toward the trailing edge of the fan. Although computations 
in the expansion fan have not been extended past | = 0.62 = £ z * in the 

fan, two prime possibilities exist for flows with values at the trailing 
edge of | > 0 . 62 . One is the shock-free flow diagramed in figure 28(e) 
in which the values of | for all the Z- characteristic lines approach 
lz* = 0.62 asymptotically. The other possibility, sketched in fig- 
ure 28(f), shows a Z-shock which must be located at £ < 0.62 to permit 
characteristic intersection since the characteristics from £ = -1 can- 
not cross this line. Until the hyperbolic differential equations are 
solved for this region, the choice of one or the other of the flows would 
be purely speculative. 

The method of procedure for solution where |^. e > 0.62 would be 

similar in many respects to that employed for determining the flow and 
shock location of a supersonic field about a cone (that is, integration 
backward toward the axis until certain conditions are satisfied on a ray 
from the origin of the conical system) . 

For the shock-free case, a possible procedure would be: 


(a) Assume a value of 



at | = 0 . 62 . 


at 


( b ) Compute 



I 


0.62; W and 


i = 0.62 obtained by integration from £ 


Z are equal to the value 

= -1. 


(c) Integrate forward / > 0 ) along ( — ) and backward /b| > 0 


1 6 t > 0 , 


6 ! 


. St < Oj 


along 


'6Z' 


to obtain values of W and Z at various values of |. 


(d) By trial and error, find the value of | at which the trailing 
edge could be located in order to satisfy simultaneously equations (ll 4 b) 
and (ll4c) with the values of W and Z obtained in (c) and also the 
value of Z obtained by integrating forward from the corresponding shock 
to the trailing edge. 
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For "the shock case the procedure would be slightly more complicated: 

(a) Assume a Z-shock location 


(b) Solve simultaneously equations (ll4b) and (like), with 

a e 

for A, b, and 0 on the positive side of l^s 

(c) , (d) Proceed as in (c) and (d) for shock-free case. 

Two points of interest in regard to the preceding outline of pro- 
cedure should be noted. First, it is not to be expected that every assumed 
slope at | = 0.62 or shock location l^s would yield an eventual solu- 
tion; in fact, one of the cases might be entirely eliminated and only a 
portion of the remaining case would be valid. Secondly, the possible 
difficulty which would be introduced if the W-characteristics either 
merged or reached their asymptotic value is not considered. By employing 
the anchor-point procedure to locate it was found that, for a shock 

tube employing air throughout, this difficulty arises at | w * = Ite = 2.25. 


THE COMPLETE SOLUTION FOR THE BOUNDARY LAYER IN 
A SHOCK TUBE FOR 7 = 1.4, R e = Roo, AND A ^ 20 


Combination, on the basis of the boundary-layer characteristics 
system, of the component solutions described in the foregoing sections 
results in the general solution for the laminar boundary layer in a shock 
tube in the entire region affected by the shock or expansion wave. The 
solution for the case of | te = 0 is illustrated in figures 29 to 52 by 

the solid lines. The dashed lines show possible smoothing or rounding 
off of the discontinuities. The trends shown by these figures are typical 
of the other values of l-^e* 

Since the solution for the flows inside the expansion fan and between 
the shock and entropy discontinuity have been discussed in detail, a brief 
discussion of the behavior only in the region | te < | 5? is presented. 

Here Z was computed analytically and W, by graphical integration. Fig- 
ures 29 and 50 show Z to be linear and continuous in this range with a 
discontinuity in slope at | 2* j whereas W is almost linear for 
l^ e < | < and then displays very strong curvature as I -* Id because 

of the adjusting process required after the entropy discontinuity. The 
value of b changes very little for $te < t < £w* 8X1(1 then becomes ver Y 

small as W becomes large (£ 1^) . 
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In the constant-pressure region between and | z *, and qg 

decrease in a manner similar to the decrease found either in the flow over 
a semi- infinite flat plate or behind a shock or negative shock as the dis- 
tance from the leading edge or shock increases. In fact, in regard to Z 
and t £ , the flow between and is identical to that behind a 

negative shock of velocity Is = ^ZO (igo ■'• n ^8* 29 where 

linear extrapolation of the curve of the' variation of Z with | inter- 
sects the Z-axis), the outer potential-flow conditions being identical 
to those in region £. 

The value of r e reaches a minimum at and increases continu- 

ously and monotonically thereafter to a theoretically infinite value at 
the shock. On the other hand, rises slightly from | z * to £ w * 

and then begins the rapid adjustment to the large negative values existing 
at the entropy discontinuity. The zero heat-transfer point appears in 
the interval between £ w * and | d - Between the entropy discontinuity 

and the shock, the value of decreases continuously and monotonically 

to -00. 

A composite diagram of the distribution of velocity and enthalpy in 
the boundary layer of the shock-tube flow is shown in figure 33 to com- 
plete the description of the entire boundary- layer flow for i te = 0. 

In parts (a) and (b) of figure 33, the profiles existing at the various 
values of | are plotted at that value of £, and in part (c) the extent 
and location of the various waves and regions are indicated. The profiles 
existing on both sides of the trailing-edge discontinuity are drawn. The 
enthalpy profiles existing on opposite sides of the entropy discontinuity 
have been illustrated previously in figure 25. The maximum boundary- layer 
thicknesses are located at the entropy discontinuity at any given time t 
and a second smaller maximum in velocity boundary- layer thickness is 
located at 

Complete distributions of t € and in the shock tube for various 

matched values of expansion- and shock- wave strength are plotted in fig- 
ures 34 and 35 . (Large-size working plots of figs. 34 and 35 are available 
on request from NACA Headquarters, Washington, D. C.) These curves, similar 
to those for |^- e = 0 discussed previously, permit the evaluation of the 

wall shear and heat transfer at any point x,t or £,t in a shock 
tube by using a single gas having 7 = 1.4 and Np r = 0.72 for shock- 

pressure ratios up to approximately 5* The dashed line, representing 
the expansion-fan solution, is used for £ ^ Ite and the solid lines, 

employed thereafter, are used for the constant-pressure regions £ and a. 
These curves may be used in conjunction with the equations of reference 7 
to predict waves generated by wall effects on an averaged basis (thick 
boundary Dayer). 
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If the method of reference 11 (thin boundary layer) is used to pre- 
dict the wave generation, it is necessary to determine the vertical veloc- 
ity Vg existing at the edge of the boundary layer. Integration of the 

continuity equation (eq. (B3)) from y = 0 to y = 5 and application 
of the boundary condition v w = 0 results in the following expression 

for vg: 



The nondimensional parameter 



is plotted as a function of 


I 


in 


figure 36 for the expansion fan in general and for particular values of 
expansion-fan strengths = -0.6, 0, and 0.5* Also shown in this 


figure is the empirical curve 


v 5 


\fl. 


1 + l< 


This curve gives 


a very close approximation to the computed results in the fan and would 
prove to be useful in the analytical integrations required if the method 
of reference 11 were to be extended to consider the expansion fan cor- 
rectly. The rapid increase in vg as | -» £4 is due to the rapid 
increase in boundary-layer thickness in this region because of the pres- 
ence of the hotter gas (low density) from a entrained in the boundary 
layer. 


APPROXIMATE SOLUTIONS FOR THE BOUNDARY LAYER IN A 
SHOCK TUBE USING ONLY ONE FLUID 


The complexity of the solutions to the reduced hyperbolic differential 
equations for the shock-tube flows is such that an approximate solution 
resulting in a large reduction in effort with only a small loss in accu- 
racy would be most desirable. Such an approximate solution is described 
in this section. For purposes of discussion, the previous numerical and 
graphical solutions to the differential equations will be termed "correct" 
or "exact." 
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Negative -Shock Approximation for the 
Expansion Fan and Region £ 

The simplest approximation available is the replacement of the expan- 
sion fan by a negative shock. This substitution was first suggested in 
reference 7 and- was also used as the basis for the computations of ref- 
erences 11 and 13 • In figures 29 to 32, which represent conditions in 
an air-air shock tube with = 0, the results obtained by the negative- 

shock assumption with both consideration and neglect of the critical char- 
acteristics axe plotted for comparison with the correct integral solution. 
(Similar results are obtained for other values of the discrepancy 

between the negative-shock and correct solution increasing or decreasing 
as increases or decreases.) 

The skin-friction parameter t £ in the region £ ^ 1^* does not 

depend upon the expansion when the critical characteristic is considered 
and is hence the same for both methods; this fact is also true for 

in the region £ > It is evident that the negative-shock distri- 

bution is a very poor approximation to the skin friction and heat transfer 
within the fan and that qualitative agreement exists only in the region 

Ste S I £ S z *- 

If it is assumed (by neglecting the critical characteristics) that 
the negative -shock solution is valid from the leading edge of the expansion 
fan all the way to the entropy discontinuity (refs. 7 , H, and 13), the 
distribution between < | !> is shown by the heavy extensions to 

the negative- shock curves. Neglect of the critical characteristics is 
tantamount to assuming that the flow generated by a negative shock which 
originated at time t = 0 is equivalent to a flow generated by a nega- 
tive shock of the same strength originating at t = -«>. 


Modified Negative -Shock Approximation 

Inside expansion fan .- In view of the failure of the simple negative- 
shock method to produce results that even qualitatively approximate the 
true expansion-fan results, a more refined method of approximating the 
expansion-fan solution with the negative shock as a basis was devised and 
applied to the entire shock-tube boundary layer. The relation between the 
true fan and the negative -shock wall-shear and heat-transfer functions for 
the vanishingly weak wave has been discussed and was found to be: 
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For I = 4te “ 1: 



Fan 



Negative 


shock 


2.000 


and 




Fan 


\(He)] 


= 2.000 


Negative shock 


If these relations are assumed to be valid not only in the limit 
t _i but also throughout the entire expansion fan, the distribution 

of 


A 


DUU cLJ_b U OIUUU^UUUU wnv. ^ jr ' _ _p 

and ^ may be approximated as follows: At a given value o 
li ™ for lot. t.ha value of t_ or q^. be twice the value 


inside the expansion fan, let the value of r e 
(at this i) obtained from the negative-shock solution for a constant 
potential flow identical to that existing at |. 

The results of this approximation along with the correct expansion- 
fan results axe plotted in figures 37 and 38 for 7 = 1.40 and Nj* - O.72. 
The approximation for r e is seen to have a negligible error for | < 0 
and an error of only 13 percent at | = 0.5 (approximately the critical 
characteristic for the negative shock). The approximation for diverges 

at lower values of |, reaches an error of 20 percent at | = -0.5, and 
increases thereafter. The reason for the relatively close lies 
t over a large range of | and the more rapid divergen 

in the fact that is much more sensitive to the boundary- layer past 

history of varying values of % at the edge of the boundary layer. In 
fact, the approximate value of % would not become zero until after 

Hi > H v (r < 1) whereas the correct q e becomes zero at | " 0-5 where 

Hi < H w (r > l). Therefore, in order to compensate for this effect, the 

approximate curve for the variation of % with £ in the range 
0 5 < | < 0.5 is replaced by a straight line originating at the approxi- 
mate value of % at | = -0.5 and passing through zero at I = 0.5- 

The equation for Hp* is 


H-l* = 21Jl2(|2 - DH* 
7+1 
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Consequently, the value of Hq* passes through a minimum at £ = 0 

imation°should bf 1 = ^ + irr fP ective <* W therefore, the linear approx- 
imation should be appropriate for other values of y. The possible error 

introduced by the linear approximation for values of Np* other than 0.72 
is not predictable . 


— region in the limit of | -> -1, the relation 




Negative shock 


O.56 


1 

2 


for ^ He was found to apply, it is now assumed that this rela- 


tion is also applicable at the trailing edge of all finite 
In addition, at the trailing edge of the fan, the matching 


expansion fans, 
solution showed 


that the approximation « 1.0 was appropriate. Thus the value of Z 


on the £-side of the trailing edge of the 
approximated as 


expansion fan may be closely 


Z(£te) , 


= |[z(Ste)] 


Negative shock, £ 


(-1 £ 6 £ 0.5) 


(147) 


where the negative shock is that which produces conditions of region t 
ehind it. Figure 39 shows that this approximation compares closely with 
the correct values. The Z-distribution in the range £ te ^ ^ | z is 

linear and may be computed from equation (102). The shear stress follows 
from equation (130) if a zero value for A is used. 

The heat-transfer parameter % for 5 te = < i z * was computed 
from equation (133a) by using the values of described previously 

and values of b(Bi) from the negative- shock solution. In this respect 
it should be noted that the value of b in the correct solution changes' 

on y s ig y m the interval |q e S £ £^* and, furthermore, the value 

of b for the corresponding negative shocks is a good approximation to 
the b-values in the region £ for £ < £ w *. 
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The distribution of for l z * s: I £ l w * was assumed to be linear 

(correct solution was nearly linear, fig. 35) a 11 *! 4he value of (^(iy ) 
was computed by the anchor-point method. For the interval ly* ^ ^ = ^d* 
an inspection of the correct solution shows the following: 



<U 


0 


at I = l w *> 





« 


-1 


at % = la, and 

(^e) ^ a 

at | = l d - 


Consequently, the following interpolation formula which gives 

=0 , (^) = and (y = (qe) a at | d was used (the 

\ d W| w * \ d| A d 

square-root term was used to make the slope a weak infinity at | d ) : 


Vl) = qe(6 w *) + < 


y id) 


a 


\ 


- y*w*) 




(148) 


The value of 



at 


is known from the positive- shock solution. 
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Summary of modified solution .- The refined approximate solution for 
the boundary layer in a shock tube when the same gas is used throughout is 
outlined as follows: 


Inside fan: 


[^e(l)] * 2 f ? e(lte)l (-1 £ I £ 0.5) 

Fan Negative shock, 

jjleU)] « 2 [q€(ite)J 

Fan Negative shock, | te =| (-1 £ £ g -0.5) 


[leCoJ - 2(0.5 - I) [qedte)] 

Fan J Negative shock, lte=-°-5 

(-0.5 ^ I ^ 0.5) 


Region £: 

( a ) £ = !te : ^(^te) 1=3 ~[z(lte)1 

2L - 1 Negative shock, | te 

(b) £te = £ = lz* : The value of Z is obtained from equation (102) 

by using | te as the base point; t £ is obtained from equation (130); 

and q^ is obtained from equation (133a) by using values of b and 
from the negative-shock solution. 

( c ) ig* - 5 = The value of Z is obtained from equation (102) 

by using | d as the base point; T e is obtained from equation (130) 

(Z and r e are also exact values); q^ is linear; and the value of q € (| w *) 
is obtained by the anchor-point method. 

(d) l w * ^ I £ ! d : The values of Z and r € are the same as those 
for condition (c) (these are exact values); and qg is obtained from 
equation (l48). 
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The results of this approximation for the shock-tube flow are com- 
pared with the correct solution for three expansion and shock strengths 
(7 = 1.4 and Npr = 0-72) in figures 40 and 4l. In region cr the 
approximate and exact solutions are identical. The agreement between the 
approximate and exact methods is seen to be very favorable for the cases 
represented. 

This approximate method is assumed to be applicable for Npr other 

than 0.72 and the approximate wall-shear and heat-transfer distributions 
have been computed for the case where the same gas exists in both e 
and a> with 7 equal to 1.4 and Npj. equal to unity. The appropriate 

shock and negative- shock results were utilized and heat- transfer results 
are shown in figure 42. The skin friction is identical to that for 
Np r = O .72 for the preceding approximation since the constant potential- 

flow momentum equation is independent of Npj- . However, in the exact 
solution Tg depends on Np^- in the fan through the pressure-gradient 

term in the momentum equation; this dependence would influence the entire 
region | g | z *. This dependence is probably small and is neglected 

herein. 


APPROXIMATE SOLUTION FOR THE BOUNDARY LAYER IN A SHOCK TUBE 
USING TWO DIFFERENT FLUIDS 


In order to obtain large shock-pressure ratios in constant-area 
shock tubes, it is common practice to use a gas with a high speed of 
sound as the driver gas in the high-pressure chamber. Hydrogen or helium 
are two such gases often employed. When various expansion-fan strengths 
for hydrogen or helium are matched to the appropriate values of shock 
strengths for air, the boundary- layer solutions of the resulting shock- 
tube flow in the regions | < 0 or 4 > Id be determined by either 

the exact or the approximate methods described previously. 

However, the region 0 ^ 4 ^ 4a presents additional difficulty if 

an exact solution is to be found. This difficulty arises because the 
boundary layer in this region contains both gases which were initially 
ahead of and behind the diaphragm, the inner part of the boundary layer 
being comprised mainly of the gas that was originally in c and vice 
versa. The relative distribution and concentration of these gases through 
the boundary layer will depend on many factors, one of which is diffusion 
normal to the wall. Instead of attempting to solve this problem exactly, 
the approximate method described in the following paragraphs was employed. 
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First, it was assumed that the wall effects of the boundary layer in 
the region 0 ^ ^ 4^* were identical to those which would have existed 

if the boundary layer contained only the gas of region This descrip- 

tion of the boundary layer becomes more erroneous as 4 increases from 
zero, since the relative concentration of the gas 0 will increase in 
the interval from | = 0 to An approximation to the magnitude of 

the errors so introduced can be obtained from the following sketch in 
which the velocity profile for A = 0 is' assumed to represent roughly 
the locus of the undiffused contact surface in the boundary layer between 
the diaphragm station and the entropy discontinuity. 



At t * (— = — = 0.404 j, the contact surface inside the boundary 

z \u x 6* / 

layer is at a value of a « 0 . 17 ; thus, for 4 < most of the boundary 

layer is comprised of gas 4* 

The heat transfer and skin friction for 0^4^ 4z* can now be 

found from forward integration along the characteristics originating 

in 4 < 0 or by the modified negative-shock approach discussed previously. 


The next assumption was that across the contact surface the momentum 


deficiencies (that is. 


r^- 


u 

~1 


dy ) of the boundary layer were 


equal, a known fact for the case of the same gas on both sides of the 
discontinuity since 0j_ = 0j. This assumption appears to be reasonable 

because the velocity boundary layer in the vicinity of the entropy dis- 
continuity is comprised almost entirely of the ga6 a. 
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Since the momentum thickness 0 as well as Z in region a is 
usually computed by using values of p w for the gas a, the necessary 

general conversion factors are summarized here. In this section, sub- 
scripts denote the reference condition, and the conversion is based on 
the equivalence of the momentum deficiencies at any | . 


Equations ( 1^9) and (150) would be used, for example, for a hydrogen- 
air shock tube to evaluate Z in terms of hydrogen in state e from a 
known Z in terms of air in state 00 . 

A linear variation is now assumed as the simplest approximation to 
the distribution for Z € and P W P W in the range ^ - ^d> 'thus 





(150) 



where 


Mid) = Mid)— 


and 





(152) 
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since T w = T e = and since the boundary- layer gas at the wall was 
assumed to be £ at and a at The shear stress may then be 

evaluated from equations ( 130 ), ( 151 ) > and ( 152 ) as 



2 

1 - lz* 

Pw( £d.) P«> 

? e U) Py( £ ) n w U) ( 

/z € Uz*) 

+ - *z* 

Pw(^Z*) 


Ze(l) l 

L 1 ' 12 * 

1 £a - *z* 

z e (!a) 1 

z e (l z *) 


(153) 


For heat-transfer computations, the region g. | g is divided 

into two areas by the line |y*, its value having been determined from 

the value of b of the appropriate negative -shock solution. The anchor- 
point method of the previous solutions is not applicable here because Z £ 

is not unchanged in the interval 1^* = I = Is when the region a is 

shrunk to zero width. However, the values of b found in the other 
solutions from the negative-shock and anchor-point methods did not differ 
greatly; therefore, the value of (which is a function of b) 

obtained by using b of the negative-shock solution should be a fair 
approximation. The value of q^ly*) is then determined from equa- 
tions (133a) and (153) with the use of the negative-shock values of b(Bq) 

Again, as in the single-gas approximation, a linear variation is used 
for q^ in the interval g | g | w *, and the interpolation formula 

(eq. ( 148) ) is applied in the interval |y* ^ I ^ Id* The value of the 

heat transfer at the entropy discontinuity to be used in equation (l48) 
is 


QeUd) 


lco(!d) 



(154) 


The summary of the modified solution for two different gases may be 
outlined as follows: 

(a) | < or | > 1^: The solution is the same as the corre- 

sponding solution for a single gas. 

(b) g £ g The value of is obtained from the negative 

shock solution; the value of r £ (|) is obtained from equation ( 153 ) J the 
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term q^l) is linear from q e (|r/*) 4° and the value of 

le(lw*) ls obtained from T € (|y*) and equation (l33 a )« 

( c ) £y* = ^ = £d : The value °t T e (|) is obtained from equa- 
tion (153) j and the value of qg(i) is obtained from equation (l48) . 

Wall-shear and heat-transfer distributions were obtained when an exact 
solution was unavailable by this approximate method for two gases for 
shock- tube flows \ising air (7 = 1.40 and Np r = 0.J2) in the region °° 

and hydrogen or helium in the region e. The properties of the hydrogen 
and helium employed, the viscosity being evaluated at room temperature, 
are indicated in the following table: 


Gas 

2k 

( Np r) € 




r 00 

( Npr l 

K 


a oo 

Helium 

1.190 

1.000 

7.236 

1.082 

2.935 

Hydrogen 

1.000 

1.000 

14.365 

0.485 

3.790 


For the hydrogen-air combination, exact solutions are available 
for | ^ 3X1(1 for I < 0 ^ since 7^ = 7 air j- The assumption that 

only gas £ comprised the boundary layer for 0 S | S permitted 

the exact solution (7 = 1.4) to be extended to for this case. 

The approximate solution was used for | ^ ^ for the hydrogen- 

air combination and also for -1 S | S for the helium-air combination. 

The results of these computations are shown in figures 43 to 46 in 
which Ite ls 1X36,1 as a cross -plotting parameter. (Large-size working 
plots of these figures are available on request from NACA Headquarters, 
Washington, D. C.) Distributions of and r e desired for a specific 

shock-pressure ratio may be obtained by interpolation from figures 43 
to 46 and by use of the data from figure 47 which gives the values of 

as a function of shock-pressure ratio for helium-air and hydrogen- 

air shock tubes. 
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For helium, | z * = when |^ e = 0.435? this value of cor- 

responds to a shock-pressure ratio — « 19 for a helium-air shock tube. 

Poo 

For hydrogen, £z* = when = 0 . 508 } this value corresponds to a 

p_. 

shock-pressure ratio — « 40 for a hydrogen-air shock tube. 


CONCLUDING REMARKS 


A solution to the unsteady laminar boundary- layer flow inside cen- 
tered expansion waves and behind both centered expansion waves and shock 
waves has been obtained in this report. 

1. The general method for obtaining these solutions may be summarized 
as follows: 

(a) The unsteady Prandtl boundary- layer equations are reduced by 
integration normal to the surface. 

(b) The velocity and temperature profiles are expanded in a power 
series. 

(c) The wall temperature is assumed to be constant. 

(d) The reduced equations are transformed into hyperbolic equations 
in a conical coordinate system because the free-stream flow is conical. 

(e) The reduced hyperbolic equations are solved in closed form for 
the flow behind a shock and by numerical integration for the flow inside 
or behind the expansion fan. Multiple solutions were found inside the 
expansion fan and the correct one was selected by physical reasoning. 

(f) The integral technique is applied at the discontinuities existing 
at the trailing edge of the expansion fan and at a contact surface so that 
the characteristic solution may proceed across these discontinuities. 

2. For the problem of the entire two-dimensional, nonstationary 
laminar boundary layer in a shock tube, the following solutions were 
obtained: 

(a) An "exact" solution for an air- air shock tube (for a ratio of 
specific heats of 1.4 and at a Prandtl number of O. 72 ) was obtained with 
the method described in the preceding paragraph. Inside weak expansion waves 


NACA TN 39^ 


97 


waves and in the flow region bounded by the shock and contact discontinuity, 
solutions also exist to the complete Prandtl boundary- layer equations. 

The integral solution agreed extremely well with these solutions. 

(b) Approximate solutions in closed form, employing zero-thickness 
expansion-fan results as a basis, for the following cases: 

Air-air shock tube; Prandtl number of 0.72 

Air-air shock tube; Prandtl number of 1.00 

Helium-air shock tube; Prandtl number of 0.72 

Hydrogen-air shock tube; Prandtl number of 0.72 

The approximate (or modified negative shock) and "exact” solutions 
for air-air shock tube (Prandtl number of 0.72) exhibited very good agree- 
ment. The simple zero-width expansion-fan solutions gave very poor 
agreement . 

3. Some of the more important features of the "exact" solution are: 

(a) The existence of zones of influence limiting the extent of for- 
ward integration along a characteristic. The solution of the present 
report was arbitrarily halted inside the expansion wave at the boundaries 
of one of these zones. 

(b) The existence of momentum and energy shocks or jumps inside the 
boundary layer. 

(c) The fact that skin friction and heat transfer inside the expansion 
wave near the leading edge was twice that for a corresponding zero-pressure - 
gradient potential (that is, Rayleigh) flow. 

4. Note should be taken of the fact that the boundary- layer char- 
acteristic system in conical coordinates derived in this paper is only 
a special case of the universal compressible -boundary- layer character- 
istics system obtainable from the integrated form of the Prandtl equa- 
tions. This universal boundary- layer characteristics system is applicable 
to unsteady compressible boundary layers in general and provides a powerful 
tool for the attack on the problem of wall effects in nonstationary flows. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , January 30, 1957 
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APPENDIX A 


EXPRESSIONS FOR VARIOUS FUNCTIONS IN A CENTERED 
EXPANSION WAVE MOVING INTO FLUID AT REST 


The wave is assumed to have the coordinates x = 0, t = 0 at its 
origin or center and to have a leading edge which advances into the region 
x < 0 for t > 0. The Riemann invariant (see, for example, ref. 24) 
prescribes the relation, valid everywhere in the wave, between and 

a-^ as follows : 


% + 


- 1-1 




(Al) 


Since a characteristic (wave) moves with a speed equal to the alge- 
braic difference of the speed of sound and fluid speed, a definite rela- 
tion exists between x and t at a point dependent on the value of u-^ 

and a^ at the point x,t. 

x = / - a 1 '\a g t (A2) 


I 


x 

a £ t 


% 


% 


(A3) 


A combination of equations (Al) and (A3) results in the expressions 
for the conical form of and a^: 


% = 


7+1 


(i + 1) 


% = 


7-1 

1 - 7 -z-l 

7+1 


(A4) 


(A5) 


The assumption of constant specific heat permits the expression of 
enthalpy in terms of speed of sound. Also, since the flow is isentropic, 
the pressure and density may be related to the speed of sound. The fol- 
lowing equivalents are then obtained: 
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(A6) 


% 







2C p T e 



2 

7 + 1 
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APPENDIX B 

DERIVATION OF REDUCED BOUNDARY -LAYER EQUATIONS FOR 
TWO-DIMENSIONAL UNSTEADY FLOW 


The Navier-Stokes equations under the restriction of the Prandtl 
boundary-layer approximation may be written as 


Ou . Ou , Ou 1 Op , 1 0 0u\ 

¥ + u ¥ + v ¥ = “ p ¥ + p ¥7 


(Bl) 


0 = - i & 
P $y 


+ ^-(pu) + ^-(pv) 

Ot Ox c y 


= 0 


(B2) 

(B3) 


Oh 

Ot 


+ u 


P Oy' 


Oh 

+ V ^ - 

Ox 

¥ 

Ou^ 

1 + 

¥ 

1 pW 


+ u + u 


Ou 

dx 


+ V ¥/ = 



(B4) 


Equations (Bl) and (B2) are the momentum equations in the x- and 
y -directions, equation (B3) is the continuity equation; and equation (b 4) 
is a form of the energy equation that is desirable when the local enthalpy 
is to be a parameter of the solution. An alternate form of the energy 
equation employing the total enthalpy as a parameter is given by 


& + uM + v aB.iA( u SB) + liL 

Ot Ox Oy P 0y\ Oy I POy 


Np r 


t 1 ' **>(1 ' u I) 


+ 1 cfc 
P Ot 


An immediate consequence of equation (B2) is 

p(x,y,t) = p(x,t) = P]_(x,t) 


(B5) 


(b6) 
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Combining equations (b 6) and. (Bl) evaluated, in the free stream yields. 


= -p,(^ * 

ck Sx 1 \St 


Un 



(bt) 


The technique of integrating each term from y - 0 to y 00 
(see, for example, ref. 25 for steady-flow integral method) is then 
applied to equation (Bl). The velocity component v is eliminated 
through substitution of the continuity equation (B3) and the following 
boundary conditions are employed: 

At y = 0 for all x,t 

u = v = 0 


(B8) 


The following momentum equation results: 


At y = 00 for all x,t 


At y = 00 for all x,t 


u = Un 


<*i = 0 

$y 


dx 





p , p ^1 

— dy + — ^ — 

~e U 1 





dy + 


1 ^1 , 
U 1 dx 



f l - JL 
( U 1 


— dy + 

Pe 


U 1 




/ du \ 

P e U l 2 Ww 


(B9) 


Equation (B9) is then rewritten by substituting 0 and 5* into equa- 
tion (B9): 
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» + -L^e + 

Sx u i Sx 


ll_ ^1 + J_ ^l\ 

b* + f (— - l)^ dy 

\ u i Sx Ul 2 st y 

Jo \ p / p e 


j__ _ ^w (S u 

u l St PgU^W. 


(BIO) 


w 


For a perfect gas with a constant specific heat, the integral term in 
equation (BIO) may be evaluated as follows: 


rffi-iW- r [iL_i\^dy = _^r (i-2L\-e_ay 
Jo \ P /p£ Jo \ h l Jp e h l J 0 { h l*J P e 


hp* * 
h l 

(Bll) 


or 


0 


Pi 

P 


1)^ dy 


H 1 “ 


ui 2 L 


Hp*a* + — + e ) 


(Blla) 


Substitution of (Bll) into (BIO) yields the momentum differential equa- 
tion with x and t as the independent variables: 


l + lh(s* - h* v4 ♦ i ||i = (BIS) 

• ^ U 1 Sx U p 2 Dt l h L j U 1 St p eUl 2 Ww 


or 


Se 


2_ l^i 0 + _i_ lii 


Sx u, Sx u 2 Dt 


6 * + 


u, 2 

-Hp*rj* + — (b* + 9 ) 


Hn 


u 


_1_ Se>5 _ 

Up St 


P e u L 


W I Su 

2 W, 


w 


(B12a) 
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The energy equation (b 4) in terms of local enthalpy is treated in 
a similar mannerj again the continuity equation is used to eliminate v 
and the following free-stream equality is substituted: 

= + = (B13) 

Dt Dt 1 Dt P L dt 

Appropriate boundary conditions for the local enthalpy are: 


At y = 0 for all x,t 


At y = 00 for all x,t 


h = tv 


h = hn 


At y = oo for all x,t 


= o 


The resulting equation for a perfect gas of constant specific heat is 


A r uA . fcMiL ay + Al !a. + .a. !h! . _k- x 

^xj o u ]_ \ ^17 \1 ^ ^1* Pi^i ^ / 




i / 1 ^i 

u il p i h i ^ 


^^!U7i 
v )J o ^ 




i ( ah i 

u lhi*\St 



+ 




— dy = 
P € 


1 ky /dh\ 

P e U l h l* c p,w vWw 



(B14) 


The last term on the right-hand side of equation (Bl4), which repre- 
sents the viscous dissipation, may be manipulated in the following manner. 
First, the Chapman -Rubes in relation (ref. l8) between viscosity and 
temperature 
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^ _ r 1 


-ft) 


1 

\l/2 


T-i + S 
T w + S 


is combined with the gas law and equation (b 6) to give 

pq 

Pw^w 


= 1.0 


(B15) 


(Bl6) 


a relation valid across the boundary layer at arbitrary x and t. 
Also since 


and 



(B17) 


da _ 1 p 
dy A p e 


(Bl8) 


d d -SL 

du _ u U 1 da _ ^1 _g_ U 1 
dy ^ da dy A P £ da 


(B19) 


Therefore, 





^w^w 

A 




2 


da 


(B20) 


d iL 

since - = 0 for a = 1.0. Equation (Bl4) is now rewritten to give 

OCX 

the energy equation in terms of the local enthalpy with independent 
variables x and t: 
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M + ( 


Sx ' 

1 — 1 

1 I 



\st i 


Shp* p Sp 


2i)d - M-L. - -i- |^)v* - 

Sx r u i\Pi h i <*t h i* / 


^Pl\ c „ ^ 1 ST* _ 1 fSh 


Up St P € u i h l* c p,wWw u l h l* 


(B21) 


Similar integrations are performed with equation (B 5 ) "by utilizing 
the following boundary conditions on total enthalpy: 


At y = 0 for all x,t 


At y = 00 for all x,t 


At y = 00 for all x,t 


H = Hw 


H = Hp 


, 0 


dy 

The following equations, which are the total enthalpy forms of the energy 
equation and are equivalent to equations ( Bl4 ) and (B2l) result. 


<L r - P,\iL ay + 1 + 1 3Hl *' 


SxJo Up l Hp* i p € 


— (l - Ej- 2 - dy + 


Up Sx Hp* Sx ) J 0 Up l Hp */ p e 


jHj* ot J 0 y H 1 *jp s Hj_* dx J 0 [ UjJp 6 
J 0 \P / P e Up St J 0 \ h 17P g 


P 1 U 1 H 1 * ^ 

1 p w ^Sh\ 
u l H l* p e ^PTjwWw 


(B22) 
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d* + fi 

dx l u ]_ dx H]_* dx 


dH,* 


±1 

1 1 

u l( 

^-4 


dp 


dH,*\ 


1 y l 1 U11 l \ * 


2 Pl dt H x * dt 


p + 


Un 


1 l u , li + __? -1 i±W* + 1 


U l' 


~ 5 ~ 1 


u^H^* \ ^ dx 


H l-^ 


2 p, dt / u^* u 2 P x dt 


1 . fPi e + I^ = 


Ul 


dt 


H n - 


1 " ~ 


1 Pw ( dH\ 


u 1 H 1 # ^e »Pr,wWv 


(B23) 


The reduced partial differential equations in final form for x 
and t as independent variables are: 


de 

dx 


JL ^ + _JL_ ^lfe* 

U 1 ^ ui 


h n * 


2 Dt 0 h n 


\ 

d ^ 

l^A 

+ 1 9 

b v el 

U-j_ St 

d Ji-\ 



e + 


de 

U 1 0 &t p u 2\^r 

M e 1 



(B24a) 


dx 


1 _ ^1 + _ 1 _ 


1 v* 


1 dp], _ 1 / 1 ^1 

U 1 dx h]_* dx Pl^l dx u l\Pl^l ^1* dt 


f dhi 


V* 


u^hj* \dt 


_1_ dPik V + _1_ \ 0 


Pl 


St 


A ^ u-i dt 


U]_ 0 dt 


1 k-v 


$ = 1 


p € u l h l* c P,wWw u l h l* u l\P e / N Pr,w\ S P u l h l* 



g 1 1_ 


(B24b) 
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cty 

Sx 


+ 


U]_ dx 


1 

Hx* ox 



1 1_ ^1 

U ]_ 2 Pi dt 
" IT 


1 

H]* 



+ 


u 1 H 1 * 


r u i 

u. fit —2: — i fii| c s; s + 

1 dx 


2 Pn at At 




Hi 


ui p l 


Un 


d* 


«iV 


X c 0 0 + I ^ t 


H n - 


U 1 


2 p-^ St A \J/- St 


t+ iff it 

y\r dt 


M*w /d H] 


v v/ p w \ 1 


P € U 1 H 1* "p r ,w Ww “lVV “pr,„ 



Equations (B24) are for h-j_, 0, and 0 as primary dependent 

variables whereas equations (B25) are for u-|_, H^, 0, and t as 
primary dependent variables. 


6* a© 

0 at 


(B25a) 


(B25b) 


108 


NACA TN 39^ 


APPENDIX C 


TRANSFORMATION OF REDUCED DIFFERENTIAE EQUATIONS 
TO CONICAL COORDINATE SYSTEM 


The differential equations obtained in appendix B are to be trans- 
formed from the x,y,t system to a conical coordinate system |,y,t where 



y = y 

t = t 


(Cl) 


This transformation is arbitrarily restricted to conical free-stream 
flows (that is, all free-stream quantities are functions of | alone) 
to reduce the complexity of the equations while still retaining the 
features that are necessary for application to flows generated by shock 
waves and centered expansion waves. 


The following derivatives are equivalent: 


i) t ‘ = 








(C2) 

(03) 


In particular, for any free-stream function f]_(x,t) = f-^(|), 


and 



Df- 


1 

t 


a e t/d£ 


(* - i)i 


dfi 


t d| 


(C4) 

(C5) 


Dt 


4 - 


(C6) 
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*V - 


It should be noted that hp* is not a free-stream function for 
lv(j:,t); thus the derivatives of hp* become 


^dhp* 
V dx 


" l^ h l “ ^ " a^t d^ 1 " ^ a^d^ 1 " d| V ) 


'dhp*' 


d ,, £ d /, . x + I dhl 

dt /x = dt^ 1 “ " t d|^ 1 “ ^ " dt t dg t dg 


(CT) 


Similar equations are applicable for Hp*. 

Application of these equations to equation (B2^a) results in the 
following relations: 


80 


uq dg 


dui a e 
+ 


U 1 


^^1/5* - V I 21 a\ + tsl Lx _ fd lx 

S 'dg l 9 hp b V 0/ up dt up d| 


0 + 


a e t &* de. _ fe|. 8* 80 _ ^w/£wV 
u x 0 dt Up 0 dg up\p £ / 


14 ) 

\ da 


w 


0 a e t 
A 0 


(C8) 


or 


de 


A - J- 6^ 

+ 0 2 < 

1 d % 

2 + (u n - 1 

g)-L| 

/&* 

_ v 

1 2 Aj 

|C 

t- 4 

CD 

Up dg 

l-l 


V 9 

h i 

b V 0/ 


v 5* 
t 3 TT 


% 


*N 8* 

5 3 "5“ 


dt q-i dg 


+ 


6* t0 d0 _ v w/ p w 


p /d JL\ 

, / Un | 


0 Up dt 


%\ p ey 


da / v A 


it 


(C9) 


In anticipation of a Pohlhausen power series for the velocity profile, 

Pet — = f(A). Derivatives of — may then be expressed as derivatives 
0 9 

of A G.S follows : 


^d ^ 

0_ 

dt 


d 61 


f-1 

dA \dt/ 


d6i 


(CIO) 
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Equations ( C9 ) and (CIO) are combined to give 

Se 2 
*6 





(1 - i- + 2e 2 * 

1 

0 , -1 “ 6 /s* h l* 1 '7* 

\ Si e / 

Uq <i| 

U]_ h L b V 0/ 
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_l e_( t SA 
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SA 
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s* t_ 

0 u L St 


^W| 

fP w\ 
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31 

V P €/ 
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S(X / A 
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(Cll) 


Substitution of equations (Cl) to (C 7 ) into equation (B2Ab) results 
in the following equation: 


S| 
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1 ^L + 


! dh-i 


u-^ di h^*\d^ 
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dpi a e 
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(C12) 


In anticipation of a power-series representation of the local 

enthalpy, the restriction is introduced that 2! = 2!(A,r,b). The fol- 

0 0 

lowing equations which are parallel to equations (CIO) express the 
V* 

derivatives of -7- : 
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The combination of equations (C12) and (C13) gives 

1 dp 


_ J. + 2 0 2 i i + jJp= . Js) 
S£ y Ui p j |U]_d? h L *\d £ Sg ) 


1 + d 


p i h 1 d| u x 


dp-, 


Pl h l d ^ 


1 L Sly Sly dhi' 
hi* \ dg ~ St ” 5 dg j 


¥1 _ 6 / 1 d Pl _ dh lV 6f; V 

0 u lhq*\P]_ dg d ^ J A 0 


S / 

_i LJT t sa 

u SA y St 


. i + i lift a: - 

1 51/ u ar \ at 


. a v* / 

s |e) + iAEft & 

Sg ) Sb \ St 


. Sb 
5 Si J 


1 V* ^ S(^ 2 _ j^w f Py\ 2 1_ 


u-, St u.]_ Ip 


el %r,w 



0. t 2±_ cfo 

V u^* ^ 


(cxh-) 


In a similar manner, the conical partial differential equations 
in terms of 0 and \J r may be obtained from equations (B25) as 


^1L -i-1 + 20 2 


—1 ®/ 




1_ ^1 
u-L d| 


2 + 


. IX + w ^ u * I □?: a\ 

% - 5 2 + Hl e H] - c □ e 




Si 


Hn - 




Bf, \ 

_i_ d e L sa _ , sa ' 

dA y St S g / 



(C15) 


V ui t 


2+2 







1 ^1 + 1 /A 

dHA 

+ i 

1 5 dp l , 1 / Hr 

. 1 5t . 

St d£ / 

1 u x d| + i^*\d5 

’ 517 

Hi 

2 Pl d5 Bfybt 


t£ + 

t 


W 



u i 2 A 


d l 2 

x cf / 


cLH 1 

2 5 d Pl 

1 c □ _ 

1 2 

1 dp l c 9 □ , 1 6 T( t & 

•l) 

u n 

“ 

I 

1 

1 o' 

CM 1 

ST 

1 

1—1 

W 

i 

1 a t 

W «2 

H 1 - X 

P x A i|r dA ^ dt 


q* 

_i_ a T /L 9n , i 

U L lo ^ St ' 6 U L 




u dc 

dc V dt 


H* 


, , i d* t at 2 _ 2v v A\A 2 i 1 H i* J j t 

Hi * St U L \pj Np r)W \Sx / w □ 


(C16) 
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APPENDIX D 


DERIVATION OF VELOCITY AND ENTHALPY BOUNDARY - LAYER 
POWER-SERIES COEFFICIENTS, FORM PARAMETERS, 

AND DISSIPATION FUNCTIONS 


POWER-SERIES COEFFICIENTS 


The coefficients for the velocity-profile power series may be 
expressed in terms of A2 (coefficient of a,2) and the highest order 
free-stream derivative. The values of A n are found by solving simul- 
taneously equations (ll) to (15) or (17) to (21) prescribing the boundary 
conditions on the profile together with the profile equation itself 
(eq. 8). This procedure yields: 

For the five-term series — • 


11 



(Dla) 



(DID) 



(Die) 




(D2a) 



(D2b) 
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A 4 = 5 + A 2 



a 5 



±A. 


Hi 


'3La 


(D2c) 


(D2d) 


A prime over a symbol represents differentiation with respect to the 
argument of the function in question. 

Local and stagnation enthalpy coefficients are obtained simply by 
replacing A n in the preceding equations by B n and C n , respectively. 

This substitution is permissible because the series expressions for local 
and stagnation enthalpy and the appropriate stream boundary conditions 
are identical in form to the corresponding velocity functions and boundary 
conditions . 


VELOCITY FORM PARAMETERS 


The displacement and momentum thicknesses are defined, respectively, 

by 



After transformation to the incompressible normal coordinate a through 
the use of equation ( 5 ) and introduction of the requirement that — = 1 

for a ^ 1, the relations (D3) become 
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Evaluation of these functions then requires substitution of the correct 
velocity-profile coefficients into the general relation for — (eq. (8)) 

and substituting the latter into equation (d 4) . Results of this calcula- 
tion for the five- and six-term velocity profiles are: 


For the five-term series 


u 



£ = 2JL 
A 315 


b* 

A 


A A 2 
9^5 9,072 


Jl _L_ 
10 " 120 


/ V 11 

9,072 



(D5b) 


For the six-term series — 

U 1 


6 * 1 
A 4 



(D6a) 


- = 0.10101 
A 


0.00063131A - 0.000031566A 2 + 0.0000526A ^ 


1 / 0=1 


0.0016835/'—) - 0.0000260 

' 1 / 0=1 

The definition of A = - 2 A 2 has been utilized in these equations. 


LOCAL ENTHALPY FORM PARAMETERS 


The quantity V* is defined as 


■5b 


(D6b) 



(07) 


V* = 
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After transformation to the incompressible 
the use of equation (6) and introduction of 

for P ^ 1, the expression becomes 


normal coordinate 3 through 

the boundary condition = 

hj_* 


1 



The profile function of 3 
to the term of the profile series 

the following expressions for 


° r (D8) 

( e l« (9)) and the coefficients appropriate 
are substituted into equation (D8) and 

result : 


For the five-term series - — : 

hf* 


v 10 60 


For the six-term series 




it 


p=i 


V* 

V 




(D9) 


(DIO) 


The coefficient Bg has been replaced by its equivalent r. 


Next to be considered is the quantity 0 which is defined by the 
relation 



(Dll) 


The local enthalpy boundary condition that = 1 for p > l must be 

h l 

applied in conjunction with the velocity boundary condition — = 1 for 

U 1 

a ^ 1. It is apparent that the relative thickness of the velocity and 
thermal boundary layers must be considered in order to evaluate the 
integral. For the case where the enthalpy thickness is equal to or 
greater than the velocity thickness, 5 h ^ 5 U , the integral (Dll) becomes: 



iL J£_ 
P e u-l 



dy + 
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When this equation is transformed to the |3-plane, it becomes 


Jo U A Vj Jp(M 


■ 5= • 


(D12) 


The limit 3(5 U ) is evaluated from the relations (5) and (6), and the 
following equation results: 


e(6u) ' l/ 0 U 4 dy ° V ■ b 


(D13) 


U 


The quantity b is equal to 


dy 

P e y 


and for 5 U ^ 5^, b = 1. 


f £ dy 

0 


When the velocity thickness is considered to be equal to or greater 
than the enthalpy thickness (b ^ l) , the integral (Dll) can be written as 


t-l. 


1 0 p « U1 

and is directly reduced to 

0 _ f 1 


V 


§ 1 - frb 


since /l - \ = 0 for 3^1. 


hj* 


$ 

Thus, the integral relations for ^ are dependent on whether the 

quantity b is less than or greater than unity. These integrals may 
be written 


for b < 1: 


i 

V 


I, it- fi-j* * 


> - a » 


(Dl4a) 


and for b ^ 1: 


$ _ f 1 

0 




( Dl4b ) 
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Application of the five- and six-term profile functions of AL 

T_ ^1 

and —— to equations (Dl4) yields 


For b 1, five-term series: 


%. = — i b + — b 2 y— b^ + -i- b5 + a (— . 

V 10 10 15 140 180 1 120 180 840 3,024 1 


b2 + ^ b5 


'h* 


1 + ]£ . M + sL\ + 


h l*/P=l \ 2,160 3,360 9,072/ \4-L*y p = i \ 40 90 56O 540 

/_u_\" / h b 2 b^ _ b5 \ + / u \ /h* \ / b 2 b^ b5 \ 

“ 45 140 " 432/ V u l/ a _i\ v h 1 *J 3= i V" 54o 560 " 1,296/ + 


1 

b 2 

+ j2 . 

3b4 

+ bL + 

M" A> 2 .. 

b? . b 4 

- b5 \ 

60 

45 

42 

280 

540 

W^iU70 

168 280 

1 , 296 / 




b 5 


^1,080 840 1,680 9,072 

For b ^ 1, five-term series: 

3 


(D15a) 


_2_ 1 

15 b 


l4o b 3 + 180 ^4 + ^ 


\90 b 84 b 2 560 b 3 1,080 b 5 ?/ 

*fe2_\" + 1 1 \ . 

V h i*/p = i\ 54o b 336 b 2 560 b 3 2,592 W 

/h* \" / 1 1 1 1 \ / u /1_1 J_ A 1 1\ 

V h l*/p=i\ *5 b l4o b 3 432 b 4/ b " 560 b 3 + 540 + 

/ u Y* / h* \" / 1 1 

\ Ul /c^l\ h l78=l ' 540 b 


1 1 


560 b 3 1,296 b 4 


4-1 + r 


11 11 

90 b " 420 b 5 


1 1 + / u_V' ( 1 1 _ 1 1 _ 1 1 \ ■> / 1 1 1 1 

i,512 ££ \ u l/ q-! \ 1 , 080 b 1,680 b 3 4,536 JT/ 1^1,080 b " 840 b 2 

1,680 h3 " 9,072 ^T) (Dl5b) 
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For 

t = 

V 


b 1, six-term series: 


jL _ jl + 0.1190477b 2 - 0.0297619b 14 - + 0.0138889b 5 - 0. 0021644b 5 + 

(-_jL_ + 0 ‘ 00 39^S 2b2 - 0.0019841b 1 *- + 0.001157^b 5 - 0. 0002164b 6 ) 


hf_V" (--r- - 0.0019841b 2 + 0.0014881b 1 *- - 0.0009259b 5 + 0.000l804b 6 + 

VVp^V 2 ^ 

.11 »«' / 

(}LS ( 1 ¥ L \ (-0.000066136b 2 + 0.000099206b 14 - - 0.000077161b 5 + 

\ Ul /c^iV h i / 3=1 V 

0.0000l804b 6 ) + - 0.0029762b 2 + 0.0009921b 14 - 0.0004960b 5 + 


0. 0000812b 6 ) + a(}^) (o.00004960b 2 - 0.00004961b 1 *- + 0.00003307b 5 - 


- 0. 0119048b 2 + 0. 0148809b 5 - 0.0089286b 14 + 


0.00000677b 6 ) + r 


0.0027777b 5 - 0.0003608b 6 + (— \ (-0.0003968b 2 + 0.0007440b 5 - 

0.0005952b 14 + 0.0002315b 5 - 0.0000361b 6 ) + A (0.0002976b 2 - 0. 0004464b 5 + 
0.0002976b 14 - 0.0000992b 5 + 0. 0000136b 6 ) j (Dl6a) 
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For b ^ 1 , six-term series: 

= 0.1190477 i - 0.0297619 -4 + 0.0138889 -i- - 0.0021644 \ + 
b b* b^ b 5 

fe) (-°* 001 984l J + 0.0014881 ^ - O.OOO9259 -Jj- + 0.0001804 

/Vi* \ 1,1 / 1 

(7— z] (0.0039682 i - 0.0019841 + 0.001157^ -T- - 0.0002164 i- 

\ h l /p = l V b b3 b4 b5/ 

" (- 0.000066136 I + 0.000099206 ^ _ 0.000077161 ^ + 

0.00001804 4r I + A (0.0059524 i - 0.0074404 4 t + 0.0044643 — - 

\P ) \ b b 2 b3 

0.0013889 nr + 0.0001804 -%■) + A(r 4 r\ (0. 00019841 4 - 0.00037202 -4 + 

V* 85/ \ h i / 3=1 V b b 2 


0.00029762 1 0.00011574 1 0.00001803 i\ 

? ' ? + pJ 


0.0059523 r* 

D 


it r / 

0.0019842 4- + 0.0009920 -i- - 0.0001623 — + -0.0000992 1 

b 3 b^ b5 W^iV b 

0.0000992 -4- - 0.0000661 — + 0.0000135 -4 1 + A (0.0002976 r - 

b 3 b 1 ^ b 3 y \ b 


0.0004464 -4 + 0.0002976 — - 0.0000992 4“ + 0.0000136 -4 


b 3 




(Dl6b) 
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TOTAL ENTEALPY FORM PARAMETERS 


The quantities □* and are defined, respectively, as 


□* = 


0 p e 


1 - Pr dy 


%*/ 


(D17a) 


♦ = / - 


0 


p e U 1 


H*_ 

H l*/ 


ay 


(DlTb) 


These definitions are transformed to the a)-plane through use of equa- 
tion (7). Because of the similarity between corresponding relations for 
the total enthalpy and local enthalpy, the equations for □*/□ and \|r/0 
are identical to those for the local enthalpy thickness equations (D9), 
(DIO), (D15), and (Dl6) when \|r, □*, □, c, and ft are substituted 
for "b > and r, respectively, and the stream boundary con- 


ditions 


'h* 


p=i 


and 


h* 

v;p=i 


are replaced by 


SL 

H 1%=1 


and 


HM 

K E 1*L=1 


DISSIPATION FUNCTION 


The dissipation function 0 is defined in appendix B and is given 


t>y 


0 = 


_ ^1_ PyPw 
2 A . 


r 1 


1 0 

da 


12 


da 


For the five- and six-term profiles, this integral is evaluated as 


For the five-term series 2A; 

^1 
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$ 




u l 2 p yPy 
Pe 2 " 


* 35 






-527 *(JL\" 

1,260 W/q^i 


+ 



(DIB) 


For the six-term series ^j-s 


0 = 



PyM-w 

A 






+ 
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APPENDIX E 

DERIVATION OF COMPLETE BOUNDARY- LAYER EQUATIONS 
FOR TWO-DIMENSIONAL UNSTEADY FLOW 

The Prandtl boundary-layer equations (eqs. (Bl) to (B4) ) may be 
rewritten with the aid of equation (b 6) in slightly different form as 


The "incompressible" normal ordinate Y is now introduced together 
with a stream function f which satisfies the equation of continuity 
(eq. (E2)). (A similar procedure was used in ref. 23.) 



(El) 


Sp + S(pu) S(pv) = Q 


(E2) 



Sh 

St 


+ u — + v 
Sx 


Sy St Sx 


ah, ^1 + 




Y = 



(E4) 




(E5) 
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The relationships of the partial derivatives in the x,y,t system 
to those in the x,Y,t system are: 




*>y 






X,t 


The stream function is then defined by 


u = 





(e6) 


V 


Pg /S|j 

p W/y,t 


+ 




+ 



Transforming the equations of motion (eqs. (El) and (E3)) to 
x,Y,t coordinates with the aid of the preceding equations, employing 
the Chapman-Rube sin viscosity-temperature relation (appendix B) , and 
considering the Prandtl number and specific heat to be constant results 
in the following equations: 


a 2 f . _ £1 afl = _ i + ^i Cw 

dYSt dY SxSY Sx ^^2 p dx 


(E7a) 


3h + Bh 
St SY Sx 


Sh = lfe + a* Wj 
Sx SY p \St SY Sx J 


p l^l 


C w 


.2 

1 d h 

N Pr SY 2 


+ 



(E7b) 
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The boundary conditions to be satisfied are then 


$(x,0,t) = 0 


df\ 


0 


/df 

\dY 


Y->°° 


^(x^t) 


(E8) 


h(x,o,t) = hy. 
h(x,oo,t) = hp(x,t) 


For the case of either the shock wave or the rarefaction wave 

propagating into air at rest in a uniform channel, the free -stream 

quantities become functions only of the parameter § = -2_. The equa- 

a £ t 

tions of motion may be written in similarity form by using this conical 
parameter. After introduction of the identities: 


$ = Jv ~ e t a e f(|,Tj) 

h = h e g(l,T)) 


(E9) 
Y 

v/vgt 



t = t 
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the transformation equations from the 
coordinate system are 





(&)*,* = (£),, t (£) t = 


x,Y,t coordinate system to the 



?(E10) 


Equations (E7a) and (E7b) are then transformed to the |,r),t coor- 
dinate system hy using equations (E9) and (E10) , the equation of state, 



ph, and so forth; the following equations result 





+ 


a + 

, 2 


+ 


5f\ d 2 f .11 d £i 


(Ell) 




l d 2 g 

N Pr ^ 2 


+ 


(r - i) 



+ 





+ 


SfW 


- 1 1_ d El/ df\ 

7 ^ d£ y S n / g 


7 


(E12) 
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with the boundary conditions 


f(!,0) = 0 



g(«,0) 

h £ 

g(4,°°) = = h x (0 

e 


> (E13) 
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APPENDIX F 


DERIVATION OF INTEGRAL RELATIONSHIPS AT A DISCONTINUITY 
IN THE OUTER FREE-STREAM FLOW 


In the outer inviscid flow, consider an arbitrary discontinuity 
which moves with a velocity V* relative to the wall. The mass flux 
through this discontinuity in an arbitrary distance L normal to the 
wall, L » 5, may be expressed in a coordinate system moving with a 

velocity V* as 

f L 

(Mass flux)(p e a £ ) _1 = / pu* dy (Fl) 

J 0 


The velocity relative to the discontinuity is u* = u - V* and the 
bars denote nondimensionalization by the reference density and velocity 
( Pg a G ) ’ Equation (Fl) may be evaluated by substitution and use of the 

definitions of S* and V* to give 


h * p L 

(Mass flux) (pe a e) 1 = ^ii “ ¥*) — V* - u^S* + (Hp " X*)£i dy 


(F2) 


The momentum flux may likewise be expressed as 


(Momentum flux)(p e a g 2 j =J jpu^ 


i* 2 + — — — \dy 


(F3) 


Pe a , 


ee 


(Momentum flux) (p e a £ 2 ) = _ 2V*J + ^u-^ -V*) V* 


I 0 L k * - *■) 


2 , pi 


P P a, 


e~e 


dy 


(F4) 


In a similar manner the energy flux into the discontinuity may be 
evaluated as 


(Energy flux) (P € a e ^) _1 ~J P^' 


L h + 


u 7 


dy 


(F5) 
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h, * 


(Energy flux) (P € a £ ^) = u-^ + ^ u-j^V*^ + —5= 

' o c- cL 


JV*^ - v*) 


2hi 


a. 


6 * + 


h L * 


hn 


Hl h l . Ul - y*y 

o 2 2 


V* - 1 u l 3 AI + 



(F 6 ) 


where 



da = 1(A) 


(F 7 ) 


If the discontinuity is assumed to have zero thickness and to contain 
no sources or sinks , the entering flux must equal that leaving since heat 
transfer and wall shear have zero length in which to influence the flow. 

If the conditions on opposite sides of the discontinuity are denoted as 
states i and j, the following equalities must apply: 

Mass : 



p i u i 




Pi 

e>i u i 2 


L 


(F9) 


t— IjOJ 
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Energy : 
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TABLE I.- RESULTS OF NUMERICAL INTEGRAL SOLUTION INSIDE EXPANSION FAN.. 
[BRANCH I; 7 = 1.4; N pr = O.72.J 
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ro 


1 

z 

w 

A 

r 

b 

dZ 

d 5 

dW 

as 

dA 

a? 

dT 

as 

db 

as 

A 

T e 

3 e 

- 1.00 

0 

0 
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0.8485 

0.1148 

0.1926 

-00 

00 

00 

0 

0 

-.98 

.002225 

.003760 
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-7.907 

.9038 

• 1077 

.1832 

-36.11 

52.34 

1.192 

.1305 

.1513 

-.96 

.004309 

.007337 

18.087 

-7.058 

.9220 

.1009 

.1748 

-23.63 

35-48 

.726 

.1829 

.2091 

-.91+ 

.006260 

.010755 

17.662 

-6.433 

• 9335 

.0943 

.1670 

-20.30 

28.67 

.492 

.2219 

.2498 

-.92 

.008081 

.014015 

17.278 

-5.894 

.9427 

.0879 

• 1593 

-17.51 

24.74 

.383 

.2538 

.2814 

-.90 

.009779 

.017134 

16.949 

-5.429 
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.0818 

.1521 

-16.17 

22.13 

.294 

.2811 

.3070 
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.011356 
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-5.006 

.9544 
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20.14 
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.3280 
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.0649 
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- 13.37 

17.56 

.143 

.3457 

.3602 

-.80 

.016560 

.030698 

15.559 

-3.591 
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(a) t=0 before diaphragm destruction 


^Expansion wave ^ — Contact Shock wave 
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(b)t=t|>0 after diaphragm destruction 



Figure 1.- Sketches showing various regions and waves existing in a 
constant-area shock tube before and after diaphragm destruction. 
Circled symbols indicate regions. 
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Figure 2.- Variation of skin-friction parameter as a function of relative velocity ratio for 

shock-initiated flows. 
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Figure 3.- Variation of heat-transfer parameter as a function of relative velocity ratio for 

shock-initiated flows. 



Figure 4 .- Variation of recovery factor as a function of relative velocity ratio for shock- 

initiated flows. 
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Figure 5.- Velocity profiles at | d for shock-initiated flow of 

u s/( u s " u o) = 2 *°* 



Figure 6.- Enthalpy profiles at 


for shock-initiated flow of 
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Figure 7.- Variation of skin-friction parameter as a function of 
shock pressure index for negative-shock initiated flow. 
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Ne g a t i v e -s hoc k pressure index,! 
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(a) 7 = 1.4. (b) 7 = I. 667 . 


Figure 8 .- Variation of heat-transfer parameter as a function of negative- 
shock pressure index for negative -shock initiated flow, 

7-1 
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Figure 9«- Variation of Z inside expansion fan. 
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Figure 10.- Variation of W inside expansion fan. 
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Figure 11.- Variation of A inside expansion fan. 
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Figure 12.- Variation of r inside expansion fan 
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Figure 15*- Variation of b inside expansion fan. 





Figure 14.- Variation of t g inside expansion fan. 
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Figure 15.- Variation of inside expansion fan. 
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Figure l6.- Comparative velocity profiles near leading edge of expansion fan. 
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Figure 1T»- Comparative enthalpy profiles near leading edge of expansion fan. 
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Figure 18.- Velocity profiles inside expansion fan. 
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Figure 19. - Enthalpy profiles inside expansion fan. Branch I solution 
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Figure 20.- Ratios of 0 and 0 at expansion-fan trailing-edge discon- 
tinuity. (Subscript i denotes inside expansion fan at trailing 
edge; subscript j denotes outside expansion fan at trailing edge.) 



Figure 21.- Ratios of 9 and q e at expansion-fan trailing-edge discontinuity. (Subscript i 

denotes inside expansion fan at trailing edge; subscript j denotes outside expansion fan at 
trailing edge.) 
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Figure 22.- Velocity profiles on opposite sides of expansion-fan trailing- 
edge discontinuity. Branch I solution. 
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Figure 23 •- Enthalpy profiles on opposite sides of expansion-fan trailing - 
edge discontinuity. Branch I solution. 
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Figure 24.- Velocity profiles near the entropy discontinuity. = 0; 

Id = 0-853. 
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Figure 25. - Local enthalpy profiles near the entropy discontinuity. 

^te = = 0 *^55. 
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Figure 27.- Boundary-layer characteristic diagram for air-air shock tube. 

Ite = °> = 1 » 62 . 
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Figure 28.- Boundary-layer 



momentum characteristic sketches for various values of 
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Figure 29.- Variation of Z with | in air-air shock tube. | te = Oj 

Np r = 0.72. 
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Figure 30.- Variation of W with | in air-air shock tube. 

Npp = 0.72* 
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Figure 31*- Variation of ? e with £ in air-air shock tube. = 0; 

N pr = O.72. 
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Figure 32.- Variation of with | in air-air shock tube. i te = 0; 

N pr = 0.72. 


Dimensionless normal ordinate, y 




Figure 33 .- Composite sketch showing boundary- layer velocity and enthalpy profiles as well as 
free-stream flow regions in air-air shock tube for = 0. Np r = 0.72. Circled symbols 

indicate region. 
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Figure Jk - Distribution of skin-friction parameter in air-air shock tubes. 7 - Npr - 0.72. 
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Figure 35*- Distribution of heat -transfer parameter in air-air shock tubes. 7 = 1.4; Np r = 0.72. 
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Figure 36.- Distribution of vertical velocity in gas accelerated by expansion wave 
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and approximate distribution of skin-friction parameter inside 
expansion fan. 


Figure 37* - Comparison of "exact" 
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Figure 38.- Comparison of " exact" and approximate distribution of heat-transfer parameter inside 

expansion fan. 


H 

--3 

H 


TN 39kk 


172 


MCA TN 39kk 



Figure 39. - Comparison of "exact" and approximate distribution of boundary- 
layer parameter Z in region £ for i % 1^* 
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Figure 4o.- Comparison of "exact" and approximate distribution of skin- 
friction parameter in an air -air shock tube, y = l.kj Np r = O. 72 . 


NA.CA TN 3944 


174 



- 1.0 -.5 0 .5 1.0 1.5 2.0 

Conical parameter 


Figure 4l.- Comparison of "exact" and approximate distribution of heat- 
transfer parameter in an air-air shock tube. 7 = 1.4$ Np r = 0.72. 
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Figure 42.- Approximate distribution of heat-transfer parameter in air- 
air shock tube, y = 1.4; Np r = 1.0. 




Figure 43.- Approximate distribution of skin-friction parameter in a helium-air shock tube 
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Figure 44.- Approximate distribution of heat-transfer parameter in a helium- 

air shock tube. 
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Figure 45.- Approximate distribution of skin-friction parameter in a hydrogen-air shock tube. 
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Figure 46.- Approximate distribution of heat -transfer parameter in a 

hydrogen-air shock tube. 
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Figure 47 .- Air shock-pressure ratio as a function of |^. e for various 

gases in region e of shock tube. 
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